DDC  FILE  COP’S  AEA064G79 


APIT/GE/EE/?8-l*V 


EFFECTS  OF  ANTENNA  ARRAYS  * 
ON  BROADBAND  SIGNALS  » ; 


r- 


is  LCY'*'^  THESIS  l 
• ■■■«••  *•  *—*•*•»•-•**•  •-  ■■ 

1 

Edwardyi.aL.-_,  

rayF. ttsaf 


^ ^ ^ ; 

[A  U\  AFIT/GE/EE/78  -l4)  j Edward  yi.aeka , Jr , 
r*  TJ  j '-TTST TTSll 


FEB  ! 1979 


itur 


. : V.  L 
■ w U U.I  — • 


( ili  ^ 


7<P 


Approved  for  public  release?  distribution  unlimited 

/ 


Oj.2  0^0 


f 9 01  SO  1 


} O w 


/• 


*^k>v  ..  ,,  „V 


F 

V 

V 


AP1T/GE/eV?8~1^ 

EFFECTS  OF  ANTENNA  ARRAYS 
ON  BROADBAND  SIGNALS 

THESIS 


Pioscnted  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  Training  Command 
in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science 


by 


Edward  Kaska,  Jr.,  3.S. 

Capt  USAF 

Graduate  Electrical  Engineering 
September  1979 

Approved  for  public  releases  distribution  unlimited. 


r.:;j 

lit'  ■' 

1 I- 


■'  1'UiH  SC6V.ii 

*•1.  W/W  wtSEC' 


QsJ 


i.'ij '.I'.:  vrn  'rTj)ij.ti<itnn;mn  iyr  -nnfirTTlnnr~~H!fi  Mill.  HQffMI 


< 

A. 


14 

>• 

4- 

K" 
i 


I 

k-' 

I 

«• 


Vi/ 

k 

I 

■f 

k 


1 


In  January  1978  I knew  nothing  about  antenna  arrays  or 
sproad  spectrum  signals  other  than  the  fact  that  I had  he&'rd 
the  terms  before.  Capt.  Stanley  R,  Robinson  steered  me  onjto 
a thesis  topic  suggested  by  Maj.  Jurgen  0.  Gobien.  The 
baedc  study  involved  the  study  of  adaptive  antenna  arrays  in 
order  to  determine  if  their  adaptive  nature  would  have  any 
adverse  effects  on  the  performance  of  spread  spectrum  com- 
munications systems.  By  February,  I made  the  commitment  to 
take  on  this  research  effort  as  my  Master's  level  thesis 
work. 

Capt.  Robinson  and  I spent  many  hours  together  in  the 
months  that  followed  studying  this  problem.  He  provided 
most  of  the  guidance  by  suggesting  many  possible  avenues  of 
research  effort.  I provided  the  legwork  required  to  travel 
down  each  of  these  paths.  After  having  performed  many  cal- 
culations, made  many  mistakes,  and  many  false  starts,  we 
eventually  organized  the  many  efforts  into  the  document  that 
follows . 

There  were  many  times  that  I thought  the  research  was 
going  to  be  fruitless.  Once  we  did  make  our  major  break- 
through it  became  a question  of  my  ability  to  write  faster 
than  the  clock  could  move.  The  clock  eventually  won  the 
race  and  the  effort  for  this  document  had  to  be  stopped.  As 
with  most  research  efforts,  the  material  included  here  is 
mostly  incomplete  and  just  hints  at  other  areas  for 
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continuing  research.  Also,  like  so  many  other  research 
projects,  the 'questions  answered  are  not  exactly  the  ques- 
tions we  started  out  to  answer. 

I wish  to  thank  Capt.  Robinson j my  advisor,  for  the 
many  hours  he  qlevoted  to  my  education  during  this  period. 

I also  want  to 'thank  my  sponsor,  Maj.  Gobien,  for  the  time 
he  devoted  to  get  me  started  on  this  study.  Thanks  go, 
also,  to  Maj.  Carl  and  Maj.  Carpinella  for  taking  the  time 
to  read  my  draft  and  advise  me  on  the  format  and  contents. 
To  these  people  I give  my  thanks  and  my  friendship. 

A special  thanks  goes  to  my  wife,  Elsie,  and  my  three 
daughters,  Lois,  Karen,  and  Christine,  who  stayed  with  me 
through  the  trying  times  and  long  nights  I spent  working 
on  the  details  and  the  writing.  To  these  four  ladies  I 
give  my  thanks  and  my  love. 

Edward  Raska,  Jr. 
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This  paper  develops  a complex , baseband  model  for  an 
adaptive  array.  The  array  is  assumed  to  havo  N isotropic 
elements.  Only  two  aspects  of  arrays  are  modeled i spatial 
propagation  delays  and  the  weighting  coefficients.  The 
array  model  is  used  to  determine  the  effects  of  arrays  on 
wideband  signals.  The  most  important  finding  is  that  the 
output  of  the  array  consists  of  the  input  signal  and  its 
time  derivatives.  Each  of  these  signal  components  is  mul- 
tiplied by  a complex  number  that  is  a function  of  the  array 
and  signal  parameters.  Properties  of  these  complex  numbers 
are  investigated.  A four  element  linear  array  is  used  as  a 
specific  example  to  illustrate  these  properties.  Two 
models  of  the  weighting  coefficients  are  analyzed  to  de- 
velope  information  about  "adaptive " effects.  These  models 
are  used  to  show  how  the  array  output  is  degraded  by  chang- 
ing coefficients. 


EFFECTS  OF  ANTENNA  ARRAYS  ON  BROADBAND  SIGNALS 


I Introduction 

Spread  spectrum  is  an  information  transmission  tech- 
nique whereby  a signal  is  modulated  in  such  a way  that  the 
transmitted  signal  has  a spectrum  that  can  be  orders  of  mag- 
nitude wider  than  the  information  spectrum.  At  the  receiver, 
incoming  signals  are  again  modulated  by  the  same  spreading 
signal.  This  results  in  two  effects*  (1)  all  noise  source 
signals  are  modulated  by  the  spreading  signal  and  their 
bandwidths  are  spread  ou+  by  the  same  ratio  that  the  signal 
was  spread  at  the  modulator,  (2)  the  modulation  of  the  infor- 
mation signal  a second  time  results  in  compressing  the  infor- 
mation back  down  to  its  original  bandwidth.  Thus,  when  the 
received  signals  are  passed  through  the  appropriate  filter, 
only  the  desired  signal  is  passed  without  loss.  Dixon 
(Ref  4)  is  a very  good  introduction  into  all  aspects  of 
current  spread  spectrum  technology. 

Antenna  arrays  are  simply  a group  of  antennas  that  are 
tied  together  electrically  to  yield  a resultant  antenna  that 
has  certain  desired  reception  properties . One  of  the  most 
important  uses  of  arrays  is  to  form  a highly  directional 
antenna  beam  that  will  receive  signals  from  basically  one 
direction.  Arrays  can  be  pointed  in  a given  direction  via 
two  methods*  mechanical  positioning  or  electronic  steering. 
The  mechanical  positioning  technique  has  been  used  quite 
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extensively  in  radar  applications.  There  are  several  draw- 
backs  to  this  technique i (1)  the  antenna  radiation  pattern 
is  fixed  in  shape  and  can  only  be  pointed i (2)  the  physical 
positioning  of  the  antenna  places  limitations  on  its  sitei 
(3)  the  rate  at  which  the  antenna  can  change  directions  is 
limited  by  the  mechanical  system  (4)  it  is  very  difficult 
to  place  a moving  antenna  on  the  exterior  surfaces  of  an 
aircraft. 

Electronically  steered  antenna  arrays  can  reduce  or 
eliminate  all  of  the  above  limitations.  However,  they  still 
have  drawbacks  of  their  own.  The  electronics  required  to 
shape  and  point  the  array  pattern  are  very  complex  and  ex- 
pensive. These  drawbacks  have  limited  their  use  to  appli- 
cations where  the  advantages  heavily  outweigh  the  cost, 
such  as  in  satellite  systems  (Ref  17) . 

Today  efforts  are  underway  to  design  and  build  equip- 
ment that  combines  these  two  techniques  into  a single  com- 
munications system.  Most  of  the  analysis  performed  on  an- 
tenna arrays  has  assumed  that  incoming  signals  are  monochro- 
matic or  quasi -monochromat ic . Little  analytic  work  has  been 
done  to  determine  the  effects  of  an  adaptive  array  on  a 
wideband  signal.  It  is  the  purpose  of  thi3  thesis  to  study 
these  effects.  Since  the  purpose  is  to  study  the  effect  of 
arrays  on  signals,  most  of  the  analysis  is  done  without 
interference  or  noise  being  considered  as  a parameter. 

Chapter  II  begins  the  analysis  by  introducing  the  var- 
ious models  and  assumptions  that  are  the  basis  for  the 
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ensuing  work.  Certain  spread  spectrum  concepts  are  discussed 
and  a complex  model  for  a spread  spectrum  signal  is  devel- 
oped. Next,  the  basic  properties  of  an  adaptive  array  are 
introduced.  Prom  this,  a general  input/output  equation  for 
an  N element  adaptive  array  is  developed.  Finally,  two  re- 
ceiver models  are  developed  for  use  in  later  analysis.  They 
are  a coherent  receiver  and  an  envelope  or  phase  incoherent 
receiver. 

Chapter  III  is  devoted  to  the  analysis  of  static  arrays. 
Thus,  the  first  thing  it  does  is  modify  the  models  so  that 
the  array  is  fixed.  Two  characteristics  of  the  spread  spec- 
trum signal  are  picked  out  and  studied  separately t its  wide- 
band nature  and  its  autocorrelation  function.  In  the  initial 
analysis  the  signal  is  simply  assumed  to  be  a nonzero  band- 
width signal.  Certain  properties  of  the  output  of  the  array 
are  then  developed  and  discussed.  Next,  the  properties  of  the 
output  correlation  function  are  discussed  and  demonstrated. 
Through  thiB  development,  the  effect  of  the  pulse  (or  code) 
rate  becomes  evident.  This  chapter  looks  at  three  specific 
examples!  a linear,  equally  spaced  array  of  four  elements 
and  two  four  element  sparse  arrays. 

Chapter  IV  then  moves  on  to  a study  of  adaptive  effects. 
Two  models  of  adaptive  arrays  are  studied.  There  are  a vari- 
ety of  algorithms  which  control  the  weighting  coefficients 
(Ref  16).  It  was  not  the  purpose  of  this  thesis  to  study 
these  algorithms  or  even  study  a particular  algorithm.  How- 
ever, it  was  desirable  to  determine  if  time  varying  weighting 


coefficients  have  any  additional  effect  on  the  incowing  sig- 
nal. As  a resultt  the  adaptive  nature  is  viewed  as  a per- 
turbation of  a desired  static  condition.  The  first  model  of 
the  perturbations  assumes  they  are  sinusoidal.  The  effects 
of  both  the  amplitude  and  the  frequency  are  then  investigated. 
The  second  model  assumes  that  the  perturbations  are  Gaussian 
processes  with  known  statistics  and  a statistical  analysis 
of  the  effects  is  presented. 

Chapter  V is  devoted  to  summarising  the  conclusions 
drawn  throughout  the  thesis.  It  then  proceeds  to  make  sev- 
eral recommendations  for  further  study  based  on  these  con- 
clusions. 
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II  System  Models 


This  chapter  introduces  a complex  baseband  model  for 
the  spread  spectrum  signal.  It  then  develops  a model  l'or 
the  adaptive  array.  Using  the  spread  spectrum  signal  as 
the  input,  an  equation  for  the  output  of  the  array  is  then 
developed.  This  equation  will  be  used  in  the  later  chapters 
as  the  basis  for  all  the  analysis.  Next,  the  models  for  the 
two  receivers  used  are  presented.  Finally,  the  chapter  ends 
with  a discussion  of  the  performance  criteria  to  be  used  in 
the  analysis. 

A.  Signal  and  Array  Equation 

It  is  the  purpose  of  this  section  to  introduce  those 
ideas  from  the  theories  of  antenna  arrays  and  spread  spec- 
trum techniques  that  are  used  as  a basis  for  the  main  thrust 
of  this  thesis.  It  is  assumed  that  the  reader  has  a back- 
ground in  antenna  theory  and  is  familiar  with  spread  spec- 
trum techniques.  Thus,  this  section  will  prove  to  be  a 
short  review  and  will  also  introduce  the  various  simplifying 
assumptions  used  to  keep  the  array  analysis  tractable. 

Figure  1 illustrates  the  basic  coordinate  system  that 
is  used.  It  is  a spherical  coordinate  axis  with  parameters 
(6,  r ) to  designate  position.  The  location  of  the  i'th 
element  of  an  antenna  array  is  designated  by  ( & , ^ » £ ) • 
The  source  of  the  signal  is  located  in  the  far  field  or 
Fraunhofer  region  at  angles  ( 0,  <p  ) . This  region  is  usually 

n 

defined  to  begin  at  a distance  of  2L  A from  the  antenna 
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Tc  * chip  period  i - chip  frequency 

T = word  period  c 
w 

Figure  2.  Pseudo-noise  waveform. 

where  L is  the  largest  antenna  dimension  and  is  the  wave- 
length of  the  carrier  frequency  (Refs  8:6  and  19*32-60). 

The  spread  spectrum  signal  used  in  this  thesis  is  com- 
posed of  two  parts:  the  information  signal  and  the  pseudo- 
noise signal  (PN  code).  The  information  signal  is  consid- 
ered to  he  a voice  signal  or  a digitally  sampled  voice 
signal,  v(t)  . Thus,  the  information  signal  can  he  con- 
sidered a low  frequency  signal  with  a bandwidth  in  the  range 
of  10  KHz  or  less. 

The  pseudo-noise  signal  can  he  thought  of  as  a train  of 
il's  that  switches  states  randomly  at  a relatively  high  rate 
such  as  5 MHz.  This  pulse  train,  p(t),  becomes  periodic  after 
some  fixed  number  of  pulses  has  passed.  Each  hit  of  the  code 
is  usually  referred  to  as  a chip;  the  entire  periodic  se- 
quence is  called  a code  word  (Figure  2).  Pulse  trains  can 
easily  he  implemented  by  properly  connecting  a recirculating 
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Figure  3.  Maximal  length  code  generator, 


register  of  N bits  in  length  to  create  a 2N  - l bit  pulse 
train  (Ref  4i80)  (Figure  3)*  This  type  of  code  sequence  is 
called  a maximal  length  code.  One  of  the  desirable  proper- 
ties of  the  code  word  is  that  it  has  low  autocorrelation  if 
the  time  shift  is  a full  chip  width  or  greater.  It  has  been 
shown  (Ref  9 »148)  that  the  codes  generated  by  the  recircu- 
lating register  of  Figure  3 have  this  property. 

As  Dixon  (Ref  4 » 14-23)  has  illustrated  the  power  spec- 
trum  of  a PN  code  always  has  a (sin  x/x)  character  where 
the  first  null  of  the  spectrum  is  at  the  clock  rate  of  the 
codei  i.e.,  the  first  null  of  a 5 MHz  code  is  located  at 
f = 5 MHz  . He  also  shows  that  about  90%  of  the  code  power 
is  located  within  the  first  null.  Thus,  the  code  bandwidth 
can  usually  be  considered  to  be  the  bandwidth  of  the  main 
lobe  of  the  power  spectrum.  Modulation  of  the  information 
signal  by  the  pseudo-noise  signal  creates  the  spread  Bpectrum 
signal.  Since  the  purpose  of  the  code  is  to  spread  the 
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signal  spectrum,  its  spectrum  is  usually  much  wider  than 
the  signal's  spectrum.  Thus,  the  bandwidth  of  the  spread 
spectrum  signal  is  very  nearly  the  bandwidth  of  the  code. 
This  modulation  can  be  modeled  in  one  of  two  ways.  The 
pulse  train  can  be  thought  of  as  a direct  amplitude  modula- 
tion of  the  signal  or  it  can  be  looked  at  as  a phase  modu- 
lation, , of  the  carrier  frequency,  where 

* 0 of  '7r. 

The  spread  spectrum  signal  is  then  translated  to  the 
carrier  frequency,  f0,  for  transmission.  During  any  numer- 
ical calculations  performed  later,  f will  be  considered  in 
the  UHF  range j about  350  MHz.  Bringing  these  elements  to- 
gether using  the  phase  shift  model  for  p(t) , the  spread 
spectrum  signal  becomes 

s(t)  * v ( t ) c,os  [a  v ft  i (1) 

This  can  be  written  as 

■ U> 

where  f?t{* } stands  for  the  real  part  of  the  complex  ex- 
pression. The  complex  envelope  of  s(t)  is 


ru)  • 


(3) 


Thus, 


S(t)  ■■  fit  {S(A)<i**'rt'*  } W 

J(t)  is  considered  the  complex  baseband  representation  of 
the  spread  spectrum  signal. 

The  array  is  also  modeled  usi.ng  a complex  baseband 
representation.  Since  our  primary  interest  is  to  model  the 
propagation  and  weighting  effects  of  the  adaptive  array,  a 
variety  of  simplifying  assumptions  are  made.  First,  all 
antenna  elements  are  assumed  to  be  isotropic  in  nature. 
Second,  since  we  are  concerned  primarily  with  small  arrays 
(few  elements) , mutual  coupling  effects  are  ignored. 

With  these  assumptions,  there  are  only  two  effects  to 
models  the  spatial  effects  and  the  adaptive  effects.  The 
spatial  effects  modeled  are  simply  the  relative  propagation 
time  differences  between  the  elements  of  the  array  and  the 
origin.  Typically,  one  element  of  the  array  is  chosen  as 
the  origin  and  the  propagation  delays  (or  advances)  are 
measured  with  respect  to  this  element.  The  propagation 
delay,  , is  the  time  that  elapses  from  the  moment  the 
signal  wavefront  reaches  the  origin.  Figure  4 shows  a typ- 
ical relationship  between  the  signal  and  two  elements,  one 
of  which  is  considered  the  origin.  The  propagation  delay 
from  the  i’th  element  to  the  origin  is  precisely  (X  - Y)/c 
where  c is  the  velocity  of  light.  However,  it  is  well  known 
that  if  the  angle,  a,  is  small  then  Y is  approximately  the 
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same  as  its  projection  onto  X.  Thus,  the  distance  X - Y 
can  be  closely  approximated  by  the  distance  r cos 
This  is  basically  the  same  assumption  used  to  identify  the 
far  field  region,  i.e.,  X and  Y are  approximately  parallel. 
This  approximation  is  correct  to  within  1%  error  if  the  angle 
a is  less  than  8®  . The  worst  case  for  this  approximation  is 
when  X = Y , because  r.  then  subtends  the  largest  angle  a. 
In  this  case,  the  ratio  Yj ^ must  be  greater  than  7*2  to 
keep  angle  a less  than  8*  . 

Using  the  above  approximation  for  the  propagation  delay, 
we  obtain 
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where  the  angle  can  be  calculated  as  follows  (see  Fig.  l)i 


cos  * cos  ©cos  0^  + sin  0 .sin  i).  cos(p  - fy) 


(6) 


The  propagation  delay,  ^ , has  the  effect  of  delaying  (or 
advancing)  the  signal  by  the  amount,  ^ . Thus,  the  signal 
at  the  i'th  element  can  be  represented  by 


s;  ( t)  ” s(  i * 

■ ' e*  / (7; 


Thus,  the  complex  envelope  of  the  delayed  signal  becomes 


ift(***<L>*JLirf***1 


(8) 


As  stated  earlier,  the  array  is  to  be  modeled  using  a complex 
baseband  representation.  Therefore,  the  adaptive  effects 
can  be  modeled  simply  by  a complex  weighting  term,  , 

used  as  a multiplier  on  the  signal  at  the  i'th  antenna 
element.  Multiplying  the  signal  by  the  weighting  coeffi- 
cient results  in  the  complex  output  of  the  i'th  element. 

^(fi,  ft  *)  - ?i  U)  A;  ( 4) 


•=*  3!  J J 
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Now,  since  mutual  coupling  effects  are  being  ignored  in  this 
analysis,  the  output  of  the  array  is  simply  the  summed  out- 
puts of  the  individual  elements.  Thus,  the  resultant  output 
of  the  array  is  represented  by 


(10) 

4*  ( * 

This,  then,  is  the  complex  baseband  model  of  the  array  for 
a spread  spectrum  signal.  ^ (£))  L and  (i)  are  all 
parameters  of  the  array.  They  each  have  an  effect  on  the 
output  signal  of  the  array.  Chapters  III  and  IV  of  this 
thesis  analyze  the  effect  of  these  parameters  on  the  ability 
of  the  array  to  pass  the  information  undistorted. 

B.  Receivers 

Eq  (10)  is  the  output  of  the  array.  This  becomes  the 
input  to  the  receiver.  One  of  the  objectives  of  this  thesis 
is  to  determine  what  effects  Eq  (10)  will  have  on  receivers. 
The  operation  of  two  receivers  is  analyzed  in  Chapters  III 
and  IV.  This  section  presents  the  models  for  those  two 
receivers. 

One  receiver  used  is  the  correlation  detector  shown  in 
Figure  5*  With  this  receiver  we  assume  a priori  knowledge 
of  both  the  transmission  frequency  and  the  pseudo-noise 
code.  However,  it  is  assumed  that  there  is  a start  time 
uncertainty,  T , in  the  code  generator  resulting  in  mismatch 
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a.  Pulse  String  Model 


c.  Complex  Baseband  Model 


Figure  5*  Correlation  Detector. 


between  the  transmitter  and  the  receiver.  Both  are  used  to 
modulate  the  output  of  the  adaptive  array  and  the  result  is 
passed  through  an  integrator.  This  system  would  be  used  if 
the  information,  v(t),  is  also  transmitted  digitally.  The 
integrator  period,  T,  is  the  period  of  the  information  code. 
Recalling  that  the  array  was  modeled  using  a complex,  base- 
band representation,  it  is  desirable  to  use  the  same  type 
model  for  the  receivers.  By  employing  the  PSK  form  for  the 
PN  code,  the  model  for  the  correlation  receiver  becomes  a 
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multiplication  by  the  complex  exponential,  * , 

followed  by  the  detection  circuitry  (Pig. 5c).  The  phase 
shift,  b,  is  used  to  shift  the  signal  phase  so  that  any 
desired  component  of  the  complex  signal  can  be  observed  at 
the  output  of  the  detector. 

The  other  receiver  examined  in  this  thesis  is  the 
envelope  detector  illustrated  in  Figure  6.  It  is  used  in 
an  analysis  to  judge  the  effect  of  not  knowing  or  even 
estimating  the  phase  of  the  carrier  on  the  output.  The 
baseband  model  of  this  system  (Fig.  6c)  is  simply  a corre- 
lation operation  followed  by  an  absolute  value  operation. 

Both  models  have  used  the  offset  X to  represent 
timing  mismatch  between  the  transmitter  and.  receiver  codes. 
If  X is  kept  as  a parameter,  then  the  receiver  output  is 
a function  of  t,  As  t is  allows i to  vary  the  output  be- 
comes the  correlation  function  of  the  array  output  and  the 
code.  The  properties  of  the  correlation  function  are  em- 
phasised iri  the  examples  of  Chapter  III. 

C.  Performance 

The  dos.igners  of  adaptive  arrays  consider  performance 
in  terms  of  parameters  like  the  antenna  directivity,  the 
main  lobe  beamwidth,  the  number  of  degrees  of  freedom,  etc. 
(Ref  10).  They  also  investigate  how  fast  the  array  adapts 
from  one  steady  state  to  another,  how  well  the  array  nulls 
point  sources  of  interference,  how  many  noise  sources  can  be 
nulled  or  how  many  degrees  of  freedom  are  used  up  to  null  a 
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specific  noise  source  or  group  of  noise  sources  (Refs  1.1  and 
12).  All  of  these  performance  characteristics  are  important 
and  need  to  be  known  for  any  given  system.  But  these  are 
incomplete  for  total  description  of  the  communication  sys- 
tem of  which  the  array  is  one  part. 

There  is  another  series  of  performance  criteria  that 
should  be  studied.  The  effect  of  the  array  on  a signal  can 
be  analysed  from  the  viewpoint  of  the  array  as  a filter.  As 
such,  the  array's  performance  can  be  discussed  in  terms  of 
signal  distortion.  The  array  affects  incoming  signals  in 
several  ways,  each  of  which  causes  distortion.  First,  the 
arrays  have  frequency  dependent  radiation  patterns  (as  do  all 
antennas).  These  patterns  cause  attenuations  and  phase 
shifts  of  the  various  frequency  components  of  incoming 
signals.  Since  the  effect  is  different  at  different  fre- 
quencies, this  will  result  in  a distortion  of  wideband  sig- 
nals such  as  spread  spectrum  signals.  If  the  distortion  is 
great  enough,  the  characteristics  of  the  pseudo-noise  code 
will  be  lost  and  the  receiver's  internally  generated  code 
will  not  be  able  to  synchronize  with  the  signal's  code. 
Without  code  lock,  the  system  cannot  recompress  the  spread 
spectrum  signal  back  down  to  its  original  bandwidth  and  the 
information  will  be  lost.  The  receivers  perform  a correla- 
tion of  the  received  signal  with  the  code.  Therefore,  in- 
stead of  looking  at  the  signal  directly,  we  will  investigate 
the  effects  of  the  array  on  the  properties  of  the  correlation 
function  of  the  code. 
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Second,  the  weighting  coefficients  are  time  varying  in 
an  adaptive  system.  They  are  used  as  multipliers  to  the 
incoming  signal  and  thus  they  become  a direct  modulation  of 
the  incoming  signal.  In  Chapter  IY , two  models  for  the 
weighting  coefficients  are  investigated  to  determine  con- 
ditions under  which  the  time  varying  properties  of  the 
weighting  coefficients  cause  additional  signal  distortion. 


« 
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Ill  taclgoansa  of  static  Arrays 


This  chapter  analyses  the  performance  of  arrays  when 
the  weighting  coefficients  are  fixed.  A variety  of  static 
cases  are  studied.  The  cases  analyzed  are  not  in  them- 
selves supposed  to  represent  actual  systems.  Th*  effects 
on  the  output  of  the  system  geometry  and  the  location  of 
the  signal  in  the  antenna  pattern  are  studied. 

A*  Fixed  Array  Equation 

The  analysis  begins  by  modifying  the  array  equation 
for  static  conditions.  Thus,  assume 

(11) 


where  A;  and  are  real  constants.  Then  Eq  (10)  becomes 

0(6,*,*,)  =£  *"J  (12) 


or  as  written  in  the  frequency  domain 

*e.*P)  (13) 

where  represents  the  Fourier  transform. 

At  this  point,  it  becomes  interesting  to  study  the 
array  properties  separately.  As  such,  the  specific  signal 
form  is  uninteresting.  Thus,  for  simplicity,  let  m(t) 


represent  v(j<)e^V^and  M(f)  represent  <fi(v(*)e,^p(  *}  Eq  (13) 
becomes 


o ( G , <p,  f)  = 


M(f£  * "J 

is  / 


(14) 


This  is  a function  of  the  variable  f.  o(Qt<p,f)  is  the  sys- 
tem output  and  M(f)  is  the  system  input.  From  linear  system 
theory,  the  ratio  of  the  output  to  the  input  is  considered 
the  system  transfer  function,  H(f).  In  this  case,  the  sum- 
mation term  becomes  the  transfer  function 


H(f)  (15) 

i«/  4 

However,  the  transfer  function  in  this  case  is  no  sim- 
ply a function  of  frequency.  It  is  also  a function  of  the 
weighting  coefficients  ( and  ^ ) , the  signal's  location 
( B and  <P  ) , and  the  array ' s geometry  (& , <P; , % ) • All  these 
other  parameters,  except  the  signal  location  parameters, 
become  fixed  for  the  case  of  the  static  array.  Thus,  the 
array’s  transfer  function  can  be  considered  as  a function  of 
9 , <P  , and  . The  dependence  of  the  transfer  function  on 
& and  <f  is  in  the  dependence  of  the  /j  on  0 and  <p  ( see 
Eq  (5)  for  details). 

At  this  point  it  becomes  necessary  to  introduce  the 
concept  of  the  phase  center  (Ref  19*72).  In  essence,  the 
phase  center  is  the  geometric  center  of  the  array.  If  the 
origin  is  located  at  the  phase  center,  then  the  sum  of  the 
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tine  delays  becomes  aero. 


N 

E *i  * 0 (16) 

In  an  arbitrarily  located  coordinate  system,  the  phase  center 
will  be  located  at  a time  delay  point  . It  then  becomes 
true  that 
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Thus, 
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Therefore, 
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With  this  concept  in  mind  return  to  Eq  (15)  and  factor  out 
the  phase  center  term.  The  result  is 


H(e,4>,r) . 
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Now  the  summation  term  oan  be  expanded  into  a Taylor  series. 
The  result  is 

t A a ' *a>  * 

i'l  ‘ 
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-*  A •*  6 ■*■••• 


(21) 


where 


•<»/ 


J22J 


and 


B « 


(23) 


It  should  be  noted  that  A and  B are  functions  of  all  the 
physical  array  parameters,  the  oarrier  frequency,  and 
(through  the  ^ ) the  angular  relation  between  the  signal 
location  and  the  array  orientation.  Because  of  the  func- 
tional dependence  on  a wide  variety  of  system  parameters, 
the  author  has  decided  to  make  the  functional  dependence 
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implicit  in  the  terms  themselves.  Later,  however,  functional 
relationship  with  the  signal's  location  will  be  made  explicit 
by  referring  to  these  terms  as  A( 6 ) and  B(  Q ) , where  & 
will  represent  the  signal's  location.  However,  the  reader 
should  recall  that  the  dependence  on  all  the  other  system 
parameters  still  exists. 

Combining  the  results  of  Eqs  (20)  and  (21)  into  Eq  (l4) 
we  obtain  the  system  output 

6 (e,i>,t)  -•  m(«) 

x (/*♦  -}  (24) 

The  first  exponential  is  recognised  as  a phase  shift  that  is 
functionally  dependent  on  the  relationship  between  the  loca- 
tion of  the  signal,  the  location  of  the  phase  center,  and 
the  location  of  the  origin  of  our  coordinate  system.  The 
second  exponential  is  a linear  phase  term  that  is  also  a 
function  of  the  same  points.  Eq  (24)  can  easily  be  trans- 
formed into  the  time  domain  by  recalling  that  a linear  phase 
shift  transforms  into  a time  delay  and  that  ( j*  trans- 
forms into  the  n'th  derivative.  Thus  Eq  (24)  becomes 

S(e,f  a)  . {.-(>  * <«)  e ♦ ••  • 
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(25) 
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It  should  be  explained  at  this  point  that  Eq  (25)  is  a 
baseband  representation  of  the  system.  Thus,  signals  UBed 
in  this  equation  must  be  written  in  their  comple-  baseband 
form.  For  example,  a monochromatic  signal,  cos  l+f.jb  , at 
the  array's  center  frequency,  fQ,  has  a baseband  representa- 
tion of  m(t)  = 1.0  (i.e.,  a d.c.  component)  in  this  equa- 

tion. Thus,  the  monochromatic  signal  has  no  derivative  terms 
and  the  output  of  the  array  is  just  the  first  term  of  the 
series. 

At  this  point  one  might  ask  why  the  phase  center  terms 
were  factored  out  of  Eq  (15)  before  performing  the  Taylor 
series  expansion.  The  purpose  is  to  separate  the  array 
geometry  effects  from  effects  due  to  selection  of  the  origin. 
The  terms  factored  out  of  the  array's  transfer  function  were 
dependent  on  the  selection  of  the  origin.  The  remaining 
terms  are  simply  a function  of  the  array  and  signal  geometry. 
As  Eq  (25)  shows,  the  origin  dependent  terms  result  in  a 
phase  shift  and  a time  delay. 

Now  Eq  (25)  provides  several  very  important  insights 
into  the  operation  of  antenna  arrays.  Foremost  is  that,  in 
general,  the  output  of  an  array  is  the  sum  of  the  signal  and 
all  of  its  time  derivatives.  The  signal  and  its  derivatives 
are  each  weighted  by  complex  numbers  before  they  are  added. 
These  complex  weightings  are  derived  from  the  signal  and 
array  geometry.  They  account  for  most  of  the  typically  ob- 
served properties  of  arrays.  Eq  (25)  was  obtained  through 
a power  series  expansion.  A similar  operation  can  be  per- 
formed on  any  system  transfer  function  providing  the 
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expansion  is  carried  out  using  inverse  powers  of  f,  also. 

As  will  be  seen  later,  the  utility  of  this  expansion  is  in 
the  fact  that  the  first  two  terms  are  sufficient  to  repre- 
sent the  system  output  of  arrays. 

The  first  term,  **.(i ♦ ) e has  usually 
been  considered  the  sole  output  of  the  array.  The  summa- 
tion term  is  commonly  called  the  array  factor  and  is  used 
to  determine  the  properties  of  arrays  (Ref  10 t Chapter  1). 
Such  things  as  the  radiation  pattern  of  the  array,  the  main 
lobe  location  and  beamwidth,  null  locations,  and  sidelobe 
properties  are  derived  from  this  term.  In  fact,  all  the 
usual  properties  of  arrays  have  been  studied  using  only  this 
first  term  of  the  series  expansion.  Prom  the  viewpoint  of 
this  Taylor  series  expansion  the  first  term  is  the  sole  out- 
put of  the  array  if  and  only  if  f c 0 j i.e.,  the  input  is 
monochromatic  at  the  array's  center  frequency,  fQ.  Nearly 
all  antenna  work  makes  this  assumption  or  assumes  the  signal 
is  close  enough  to  a monochromatic  signal  in  nature  that  the 
assumption  yields  valid  results. 

The  derivation  of  Eq  (25)  shows  that  the  time  deriva- 
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extremely  small  compared  to  the  first  term  and  are  not  ob- 
servable as  a result.  The  basic  reason  for  this  has  to  do 
with  the  ( ^ )~  factors  found  in  the  other  terms  of  the 

4 C 

series.  The  ^ and  ^ of  most  systems  are  on  the  order  of 
10"6  seconds  or  less.  By  attenuating  the  output  of  the 
first  derivative  term  by  this  amount  it  is  virtually  un- 
detectable in  comparison  to  the  first  term  of  the  series 
expansion  under  most  operating  conditions.  However,  in  the 
neighborhood  of  nulls  the  first  term  of  the  series  expan- 
sion becomes  almost  zero.  Thus,  in  the  neighborhood  of 
nulls,  the  derivative  term  should  be  of  comparable  magnitude 
with  the  signal  term  and  its  existence  and  effects  should  be 
observable . 

The  second  reason  for  not  seeing  the higher  order  terms 

is  related  strongly  to  the  first  reason.  No  present  day 
system  operates  with  the  signal  located  at  a null  of  the 
system.  In  fact,  most  systems  attempt  to  steer  the  array 
at  the  signal.  Recall  that  , when  £ = 90* 

the  are  identically  zero  and  there  are  no  propagation 
delays  (the  same  holds  true  for  4*  ).  This  location  is  the 
natural  main  beam  axis  of  the  array.  Replacing  the  4,-  and 
in  Eq  (25)  by  zero  results  in  the  loss  of  all  terms  of 
the  series  except  the  first  term,  the  term  containing  the 
signal  m(t).  This  means  that  none  of  tfie  signal's  deriva- 
tives are  passed  through  the  array.  The  output  is  simply 
an  amplified  and  phase  shifted  version  of  the  input  signal. 
Thus,  there  is  no  distortion.  This  is  why  the  main  beam  is 
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an  ideal  location  for  the  signal  as  has  been  known  from  all 
previous  work  in  thi3  field.  Thus,  most  operating  systems 
have  not  observed  derivative  terms  at  the  output. 

There  is  another  observation  that  can  be  made  about 
Bq  (25)  at  this  time.  In  general,  the  signal  and  all  its 
derivatives  will  each  be  multiplied  by  different  complex 
numbers.  Each  of  these  complex  numbers  will  change  the  gain 
and  phase  of  the  signal  and  its  derivatives.  Thus,  all  the 
terms  will  have  different  phase  angles.  This  means  that  if 
a phase  tracking  system  could  lock  onto  the  phase  of  the 
signal  then  only  the  projection  of  the  time  derivatives  onto 
this  phase  axis  would  be  observed  at  the  output.  In  fact, 
if  the  signal  and  its  derivatives  are  shifted  with  respect 
to  each  other  by  a full  90*  there  would  be  no  time  deriva- 
tives observed  at  the  output.  This  property  will  be  illus- 
trated in  the  example  in  the  next  section. 

There  is  one  final  concept  about  Eq  (25)  to  be  studied 
here.  It  was  shown  above  that  the  derivative  term  is  multi- 
plied by  a complex  number  that  is  usually  many  orders  of 
magnitude  smaller  than  the  complex  number  that  is  combined 
with  the  signal.  Thus,  at  most  locations  the  signal  term  is 
dominant.  However,  this  implies  that  the  signal  and  its 
derivative  are  comparable  in  magnitude.  The  true  determina- 
tion of  the  output  is  the  product  of  the  signal  and  its  gain 
and  the  product  of  the  derivative  and  its  gain.  If  the 
derivative  term  is  many  orders  of  magnitude  larger  than  the 
signal  term,  this  could  offset  the  effect  of  the  complex 
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constants  to  the  extent  that  both  the  signal  term  and  its 
derivative  term  would  be  observable  in  areas  other  than  at 
the  nulls  of  the  array.  In  fact,  if  the  derivative  is 
large  enough  it  could,  in  itself,  become  the  dominant  out- 
put everywhere  except  at  the  main  beam  axis. 

Papoulis  (Ref  15»178-182)  develops  a bound  on  signal 
derivatives  based  on  the  signal's  energy  and  bandwidth.  The 
energy  of  the  bandlimited  signal,  x(t),  is  defined  as 

E-  / I XU)I*  oU  ~ J IX(P)f  df  (26) 

'i 

where  X(f)  is  the  Fourier  transform  of  the  signal  x(t)  and 

f is  the  frequency  limit  cn  the  transform)  i.e.,  X(f)  = 0 
c 

for  If I 2.  fc  . Using  this,  linear  system  theory,  and 
the  Schwarz  inequality,  Papoulis  shows  that 

* $£>Jfc3  (27) 

The  equality  holds  at  time  t = tQ  only  if  X(f } = 
where  k is  a constant. 

This  bound  is  very  useful  in  that  it  shows  two  methods 
by  which  the  derivative  term  can  be  kept  below  certain 
desired  levels.  The  first  is  to  keep  the  signal  energy  low 
and  the  second  is  to  keep  the  bandwidth  of  the  signal  low. 
However,  both  of  these  conditions  counteract  other  desirable 
features  such  as  the  desire  to  have  E large  for  high  signal 
to  noise  ratios  and  to  have  wideband  (fQ  large)  signal 
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transmission  (spread  spectrum,  for  instance)  for  anti- 
jamming purposes. 

Eq  (2?)  is  not  an  equality)  it  is  simply  a bound  on 
the  derivative.  It  leaves  the  impression  that  the  deriva- 
tive increases  with  increases  in  the  energy  or  the  band- 
width, but  it  doesn't  prove  this.  What  Eq  (25)  has  shown 
is  that  the  derivative  term  is  passed  by  the  array.  Pa- 
poulis  has  defined  two  parameters  that  can  be  controlled  to 
keep  the  derivative  bounded  and  thus  keep  the  derivative’s 
effects  under  some  arbitrary  level.  However,  research  in 
this  area  could  determine  signal  sets  that  have  much  smaller 
derivatives  than  this  bound  indicates.  This  research  will 

I 

become  necessary  if  the  trend  toward  using  both  antenna 
arrays  and  wideband  signals  continues. 

The  analysis  to  this  point  has  shown  the  effect  the 
array  has  on  the  incoming  signal.  As  Eq  (25)  shows,  the 
output  of  the  array  is,  in  general,  somewhat  different  than 
the  signal  that  entered.  The  question  arises,  then,  what 
effect  does  this  output  have  on  the  system's  receiver. 

Eq  (25)  is  too  general  to  give  specific  results,  but  certain 
characteristics  can  be  discussed. 

First,  we  look  at  the  coherent  receiver  of  Figure  5* 

To  operate  properly,  this  receiver  must  measure  the  phase 
of  the  incoming  signal  and  set  the  local  oscillator  phase 
to  match  this  phase.  Then  the  receiver  will  be  "locked" 
onto  the  signal.  Eq  (25)  reveals  a potential  problem  with 
this  scheme.  There  are  a variety  of  signals  coming  through; 
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the  signal  and  all  ita  derivatives.  They  are  all  at  differ- 
ent phases.  The  overall  effective  phase  is  a weighted  sum 
of  all  these  phases.  A phase  tracking  system  will  probably 
lock  onto  the  overall  phase  and  thus  will  not  be  precisely 
aligned  with  the  desired  signal.  When  the  derivative  terms 
are  small  compared  to  the  signal  term,  the  composite  phase 
will  be  almost  identical  to  the  signal  phase  and  the  coher- 
ent receiver  should  function  normally.  However,  near  the 
null  points,  the  signal  term  is  highly  attenuated  and  the 
phase  will  become  more  controlled  by  the  derivative  term. 

If  the  phase  shifts  away  from  the  signal,  the  loss  of  coher- 
ence will  obviously  degrade  the  performance  of  the  coherent 
receivers.  The  example  of  Section  III.  C.  will  explore 
the  output  phase  in  detail  and  show,  for  a specific  array 
and  signal,  just  how  drastically  the  phase  can  change  near 
the  nulls. 

Although  it  is  difficult  to  determine  how  Eq  (25)  will 
affect  the  operation  of  coherent  receivers,  it  is  very  easy 
to  see  how  it  will  affect  the  envelope  detector  of  Figure  6. 
The  envelope  detector  will  just  be  the  magnitude  of  the 
overall  signal.  Thus,  if  any  of  the  terms  is  significantly 
larger  than  the  other  terms,  the  output  of  the  envelope  de- 
tector will  basically  be  that  term.  Thus,  qualitatively  we 
realise  that  in  most  areas,  the  output  will  be  the  desired 
signal.  Close  to  the  nulls  the  first  derivative  term  will 
become  dominant  and  the  output  will  appear  to  be  the  first 
derivative.  At  some  angular  distance  from  the  null,  the 
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signal  term  and  the  first  derivative  term  will  be  approx- 
imately the  same  and  the  output  will  become  their  vector 
sum.  The  example  of  Section  III.  C.  will  also  examine  this 
situation  in  more  detail. 

Although  Eq  (25)  contains  the  necessary  information 
needed  to  draw  all  the  above  conclusions  it  is  still  useful 
to  examine  several  specific  cases  that  make  the  results 
more  quantitative.  The  specific  examples  will  illustrate 
more  clearly  just  how  far  off  the  nulls  that  the  derivative 
term  is  still  significant.  The  examples  will  also  discuss 
more  quantitatively  the  effect  of  the  relationship  between 
the  signal  and  its  first  derivative. 

B.  Linear,  Equally  Spaced  Elements 

Consider  the  case  of  a linear  array  with  N elements  all 
spaced  equally  as  shown  in  Figure  7.  Also,  assume  there  are 
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no  element  weights)  i.e.i  4.  * 1 » ■ 0 . With  the 

•pacing  and  geometry  shown  tha  time  delays  baooma 

• (a  mO  i c#$  • • It  ahould  ba  noted  that  thi«  geometry 

has  removed  dependence  on  the  spatial  parameter , f>  . Thus, 
f will  be  dropped  from  the  equations  hereafter.  After 
substitution  Bq  (14)  becomes 


3(0,0 


(28) 


Poliowing  the  example  of  the  general  case,  we  would 
expand  Eq  (28)  in  a.  Taylor  series  at  this  point.  This  would 
result  in  an  infinite  series  of  summations  of  N elements 
each.  However,  the  summation  in  Eq  (28)  can  be  performed 
directly  resulting  in 


o(6,-r)  • 


tin 

c 


* irdfcta  6 1 Siftf  ? (-P+-P.  )dcos q)  /2( 

c ”*  «n[  J (f  *f9)dlcos  o] 

We  now  observe  that  M(f)  is  multiplied  by  two  exponen 
tials  and  a trigonometric  ratio.  The  exponentials  are  the 
phaBe-center  terms  obtained  in  Eq  (20)  where 


i _ AM  dec*  6 

J.  c 


(30) 
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The  first  exponential  is  independent  of  f and  is  simply  a 
complex  constant  when  the  array  parameters  are  defined.  As 
such  it  just  creates  a constant  phase  shift  of  the  array 
output  that  can  be  measured  and  compensated  for  by  a phase 
tracking  system.  The  second  exponential  is  a linear  func- 
tion of  the  frequency  f.  As  such  it  causes  a simple  time 
delay  of  the  signal  output. 

Now  the  trigonometric  term  in  Eq  (29)  can  be  analyzed 
using  a Taylor  series  expansion.  The  result  is 


»n[  ^r(-f  * fj  d c«.i  G] 
s*n[  J (f  ♦ fjrfcose] 


A * B ♦ * • * • 


(31) 


- «■»(»  r . ... 


where 


T a "n’$mClCO&  6 

' c 


(33) 


where  the  expansion  has  only  been  carried  out  for  the  first 
two  terms.  Thus, 


A , sin  (Nf) 


(3*0 


and 

m sin(N  J)o»s(J)  J 

s,na  ( J) 
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Obviously,  both  A and  B are  functions  dependent  on  the  array 
parameters  and  the  signal  location,  Q . However,  -t  .ey  sure 
shown  as  simple  constants  in  this  frequency  domain  analysis. 

A is  the  familiar  radiation  pattern  for  linear,  equally 
spaced  arrays  and  is  used  to  obtain  all  the  monochromatic 
and  "quasi -monochromatic”  results  that  have  been  published. 
Inserting  these  expressions  into  Eq  (29)  and  inverse  trans- 
forming yields 

*■)  * nn(/  + lc)  A 

* ev[  jJrr.t J g * •••  06) 

The  value  of  this  expression  is  in  the  following  obser- 
vations. First,,  m(t  + t_)  and  its  first  derivative  are 

c 

multiplied  by  the  same  complex  constant.  Thus,  this  term 
shifts  both  terms  by  the  same  phase.  Second,  the  factor  A 
is  purely  real,  while  the  factor  B,  although  complicated,  is 
purely  imaginary.  Thus,  these  two  terms  will  force  the 
signal  m(t  + tc),  and  its  first  derivative  to  be  separated 
in  phase  by  90* . 

This  is  a very  interesting  result  because  it  reveals  a 
possible  technique  for  improving  the  array  operation.  If  a 
phase  tracking  loop  can  lock  onto  the  phase  of  the  signal 
term,  then  the  phase  can  be  shifted  so  that  only  the  signal 
will  be  along  the  imaginary  axis  and  will  not  be  observed  at 
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the  output  of  the  correlator. 


The  linear  array  results  obtained  from  Eq  (36)  are 
useful,  but  quantitative  results  still  haven’t  been  inves- 
tigated to  determine  the  extent  of  these  effects.  Thus,  the 
next  section  derives  results  for  a specific  case  of  Eq  (29). 
C.  Four  Element,  Linear  Array 

The  specific  example  to  be  studied  here  will  have  the 
following  parameters! 

N « M- 

d a 

f * 350  MH* 

Substituting  the  first  two  parameters  into  Eq  (29)  yields 


ote.n  SjfoifSyf-  (37) 


This  equation  can  be  analyzed  in  exactly  the  same  way  the 
linear  array  equation  was  analyzed  in  the  last  section. 
First,  it  is  easy  to  recognize  the  first  exponential  as  a 
constant  phase  shift  that  is  a function  of  the  signal's 
location,  Q . Second,  the  other  exponential  is  the  time 
delay  term  that  is  also  a function  of  $ . It  should  be 
noted  that  both  these  terms  reduce  to  1 for  & = 90° . This 
is  the  location  of  the  main  beam  axis  for  this  array.  Thus, 
for  signals  located  on  the  main  beam  axis  there  is  no  phase 
shift  or  time  delay  caused  by  the  array*  Finally,  it  should 
be  realized  that  these  terms  do  not  distort  the  signal i they 
only  shift  it  in  time  and  phase.  However,  the  system 
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dt signer  must  compensate  for  these  terms  in  hie  design  (this 
may  not  be  trivial  for  multiple  signals  and  multiple  looa~ 
tions),  Thus,  for  this  example,  only  the  ratio  of  the  sine 
functions  is  left  to  cause  any  signal  distortion. 

In  the  last  section,  this  term  was  expanded  around 
f = 0 and  only  two  terms  in  the  expansion  were  analysed. 

No  analysis  was  made  to  indicate  when  the  third  and  higher 
order  terms  can  be  ignored.  The  usual  criteria  for  ignoring 
the  higher  order  terms  is  that  the  function  that  has  been 
expanded  is  almost  linear  in  the  region  of  interest  around 
the  expansion  point.  In  this  case,  the  region  of  interest 
is  the  bandwidth  of  the  signal. 

The  (sin  Nx/sin  x)  term  of  this  example  can  be  expanded 
exactly  like  the  similar  term  of  the  last  section.  This 
will,  in  fact,  be  done  shortly.  First,  however,  a slightly 
different  approach  will  yield  more  insight  into  when  the 
higher  order  terms  can  be  ignored. 

Figure  8 is  a plot  of  the  sinusoidal  ratio  for  several 
values  of  g . This  figure  shows  the  effect  on  frequency 
components  as  far  as  lOOJf  away  from  the  carrier  frequency, 
fQ.  For  moderately  wideband  signals  in  the  UHF  range,  the 
bardwldth-to-carrier  ratio  would  be  less  than  10#  and  in 
most  cases  more  like  or  2jC.  Thus,  the  region  of  interest 
is  confined  close  to  the  vertical  axis  in  Figure  8.  When 
f/fQ  is  small,  Figure  8 shows  that  the  curve  is  very  nearly 
linear. 

Assume  for  a moment  that  in  the  bandwidth  of  interest 
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Figure  8.  sin  (4x(e))/sir.  (x(e))  for  selected  0 


the  curve  is  precisely  linear.  Then  the  last  term  in 
Eq  (37)  could  he  replaced  hy  a term  of  the  form  A+  r . 
where  A and  3 are  real  constants  that  are  evaluated  for  each 
value  of  Q.  Letting  a , Eq  (37)  then  becomes 

3(e,r>  ■ e^rN'  (A  * ^ :~&t) 


Ae>$“*e  -.^>«6.>SwM(0e< 


3wfi% 


(38) 


Now  the  above  can  be  transformed  back  into  the  time  domain. 

S(e.y)  = 09 

The  first  term  is  simply  the  signal  shifted  in  time  and 
moved  in  quadrature  by  the  exponential  term.  The  second 
term  is  the  time  derivative  of  the  time  shifted  signal.  It 
is  also  shifted  in  phase  by  a constant  amount.  Note  that 
this  phase  shift  is  90°  less  than  the  other  phase  shift. 

If  a receiver  could  track  one  of  these  components,  it 
would  not  see  the  other  component.  If  a coherent  receiver 
locks  onto  the  phase  of  the  first  term,  the  output  will  be 
the  desired  signal.  If  it  locks  onto  the  phase  of  the  sec- 
ond term,  the  output  will  be  the  derivative  term.  If  it 
locks  onto  the  composite  phase  of  the  output,  then  the  out- 
put will  be  determined  by  the  relative  magnitude  of  the  var 
ious  terms.  As  will  be  shown  shortly,  the  signal  term  is 


dominant  in  most  regions.  Thus,  the  phase  will  be  close  to 
the  signal  phase.  Near  nulls,  however,  the  dominant  term  is 
the  derivative  and  the  phase  will  shift  toward  the  deriva- 
tive's phase  axis. 

An  envelope  detector  (as  shown  in  Pig.  6)  would  pass 
the  magnitude  to  the  receiver.  Thus,  the  output  of  the  en- 
velope detector  would  be 

«U>  .[{*»(*.*.))*♦  (40) 

after  the  quadrature  detection.  The  first  term  is  the  de- 
sired signal.  The  other  term  can  be  considered  a noise  term. 
As  long  as  this  term  is  small  compared  to  the  signal,  the 
receiver  will  be  able  to  operate.  This  can  occur  in  several 
ways i B can  be  very  small  compared  to  A or  the  derivative 
can  be  small  compared  to  the  signal. 

The  results  of  Eqs  (39)  and  (40)  were  obtained  based  on 
the  assumption  that  the  trigonometric  term  in  Eq  (37)  can 
validly  be  modeled  as  a linear  function  of  frequency.  As 
stated  earlier,  the  curves  of  Figure  8 indicate  that  this 
model  is  valid  for  f/fQ  < .1 

We  can  now  evaluate  the  constants  of  the  linear  model 
through  the  use  of  the  Taylor  series  expansion.  Substitu- 
ting the  array  parameter  into  Eq  (32)  yields 

ainfe*  ( I ♦ £ ) Cos  0 ] _ stn(3*Co&Q) 

S)n[j  (|  + £)  co60]  ‘ Sin(Sco.e) 

+ issai  rw?co .^&o,(^(>o3e)-lo^?.^e)s^X^.-0lsQLlf  (4i ) 

* L Sin  *(  J ew  6)  J 
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Thus, 


/\  * .sin  cos  a) 
sin(Jcosd) 


(42) 


and 


6 - 


4sin(  Ico*  (^7r  rn%G)  ~ Cti*.  (H60*  ?>in  (Ztt Cai$) 

0 '/'f.  Si  (V*  ( J Co  s 0 ) 


(43) 


The  term  A is  the  usual  array  factor  term  used  for 
antenna  analysis.  By  plotting  j A | in  a polar  plot,  we 
obtain  the  antenna  radiation  pattern  (Fig.  9)*  The  main 
lobe,  side  lobes,  and  nulls  are  defined  according  to  this 
term.  As  can  be  seen,  this  array  has  nulls  at  0°  and  60° . 

By  plotting  |b|  in  a similar  way  we  have  what  can  be 
called  the  "time  derivative  radiation  pattern"  (Fig.  10). 

By  observing  these  two  plots,  a very  interesting  phe- 
nomenon is  observed.  The  location  of  the  main  lobe  and 
sidelobe  peaks  in  the  antenna  radiation  pattern  is  the  loca- 
tion of  the  nulls  in  the  time  derivative  radiation  pattern. 
Furthermore,  it  is  seen  that  the  opposite  occurs.  The  peaks 
in  the  time  derivative  radiation  pattern  are  the  same  as  the 
nulls  for  the  antenna  radiation  pattern.  Later  examples 
(Section  III.  D.)  seems  to  indicate  that  the  reciprocal 
relationship  between  these  two  parameters  holds  in  general. 
No  effort  has  been  expended  to  determine  under  which  condi- 
tions this  phenomenon  exists  or  whether  it  has  any  practical 
uses.  It  is  just  pointed  out  that  the  phenomenon  exists  in 
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Figure  9*  Four  Element  Equally  Spaced  Array 
Plot  of  |A(0)|  . 


Ma*A  Value  = 2.8  x 10 


Figure  10.  Four  Element  Equally  Spaced  Array 
Plot  of  (B(0) } . 
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the  examples  used  in  this  thesis  and  the  author  suggests 
that  it  he  given  a closer  examination. 

It  ie  easy  to  show  that  J A | and  | B | are  symmetrical 
about  both  9 = 0*  and  9 = 90°  • Thus,  data  about 
these  numbers  in  the  firBt  quadrant  are  sufficient  to  de- 
scribe their  values  everywhere.  Table  I tabulates  these  two 
numbers  in  the  first  quadrant.  As  can  be  observed,  B is 
very  small  compared  to  A except  in  the  neighborhood  of  the 
two  null  points  at  0°  and  60° . Tables  II  and  III  tabulate 
these  numbers  in  more  detail  around  the  null  points.  These 
latter  tables  show  that  A decreases  and  becomes  comparable 
to  B near  the  two  nulls. 

The  values  in  the  tables  were  calculated  using  a com- 
puter to  evaluate  Eqs  (42)  and  (43).  Thus,  roundoff  errors 
are  a concern  when  the  numbers  become  small.  This  effect  is 
evident  at  the  two  null  points  0*  and  60° . By  simple  sub- 
stitution into  Eq  (42) , it  is  readily  seen  that  A is  zero  at 
these  locations.  Also,  B should  be  zero  at  9 = 90* 

The  tables  do  not  show  zero  results  at  these  points.  The 
calculations  were  run  a second  time  using  double  precision 
numbers.  The  results  at  these  three  points  became  closer  to 
the  correct  values  indicating  that  the  errors  are  due  to 
roundoff  errors  that  accumulate  during  the  calculations. 

What  is  important,  however,  is  that  none  of  the  other  entries 
changed  when  the  double  precision  calculations  were  carried 
out.  This  indicates  that  there  is  no  roundoff  error  (to 
seven  significant  figures)  in  any  of  the  other  entries. 
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These  tables  illustrate,  for  this  example,  the  state- 
ment in  Section  III.  A.  that  the  multiplier  of  the  deriva- 
tive term  is  usually  orders  of  magnitude  smaller  than  the 
multiplier  of  the  signal  term.  This  table  also  indicates 
another  valid  reason  why  the  derivative  term  is  not  normally 

observed.  As  can  be  seen  from  the  tables,  B is  on  the  order 

-8 

of  10  . Thus,  the  derivative  term  ia  attenuated  by  a volt- 

age ratio  of  80  d£  or  more  ax  all  signal  locations.  This 
amount  of  attenuation  will  normally  push  the  derivative  term 
down  into  the  system's  noise  region. 

At  this  point,  then,  it  appears  ludicrous  to  study  such 
an  insignificant  signal  further.  However,  there  are  two 
things  that  will  counter  this  actenuation  effect.  The  first 
reason  is  the  concept  of  spread  spectrum  itself.  By  using  a 
PN  coded  signal,  spread  spectrum  techniques  typically  give 
the  system  an  extra  20  - 30  dB  of  SNR  capabilities  with  a 
practical  upper  limit  of  about  70  dB  (Ref  4i7  and  24). 

Thus,  spread  spectrum  systems  have  the  capability  to  detect 
low  level  signals  if  they  know  what  they're  looking  for. 

The  analysis  of  this  thesis  indicates  what  the  form  of  that 
low  level  signal  happens  to  be. 

The  second  reason  comes  from  the  increasing  demand  for 
higher  data  transmission  rates.  For  certain  signal  sets, 
the  higher  the  data  rate  the  higher  the  time  derivatives  for 
constant  amplitude  signals.  Thus,  for  high  data  or  code  * 
rates  the  derivative  could  itself  be  orders  of  magnitude; 
higher  than  the  signal.  This  would  offset  the  80  dB  loss 
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considerably.  These  two  reasons  can  combine  to  make  the 
derivative  term  a signal  that  can  be  observed  and  can  have 
an  effect  on  the  system  output. 

The  fact  that  signals  and  their  derivatives  can  be 
orders  of  magnitude  different  leads  to  the  problem  of  de- 
termining at  which  locations  around  the  array  that  the 
derivative  term  has  a noticeable  effect  on  the  array  output. 
Tables  I,  lit  and  III  are  not  really  in  a format  that 
quickly  show  this. 

Typically  a designer  would  indicate  that  the  effect 
of  the  derivative  term  is  ignorable  if  the  signal  output- 
to-derivative  output  ratio  is  greater  than  Borne  value,  7£  . 

2*3^1 2 V 

° cm 

Rearranging  and  taking  logarithms  yields 

* s ^0®)  + t (lt5) 

Now,  the  designer  can  use  this  as  a criterion  to  determine 
those  regions  where  the  array  will  operate  satisfactorily. 
The  ratio,  log  A/B,  for  the  array  example  of  this  section 
has  been  plotted  in  Figure  11. 

To  see  how  to  use  this  figure,  a signal  and  its  deriva- 
tive are  required.  One  of  the  purposes  of  this  thesis  was 
to  determine  the  effects  of  arrays  on  spread  spectrum  sig- 
nals. This  can  be  partially  accomplished  by  determining 
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the  array  effects  on  just  the  PN  code.  We  recall  that  the 
Incoming  signal  is  correlated  with  the  code  in  order  to 
despread  the  signal.  Thus,  if  the  code  is  the  sole  input, 
we  are  ultimately  interested  in  the  array  effects  on  the 
code's  correlation  function.  It  is  true  that  the  correla- 
tion cf  a PN  code  word  with  itself  is  a triangular  waveform 
of  height  NT  and  width  2T,  where  T is  the  chip  period  and 
N is  the  code  length  (Ref  4i64-67).  Therefore,  for  the  re- 
mainder of  this  example  the  signal  will  be  considered  to  be 
the  triangular  waveform  of  Figure  12. 

Using  the  specific  values  in  Figure  12,  the  ratio 
has  the  value  5 x 10-*  near  £>  = 0 . Also,  assume  that  an 

acceptable  value  for  Tj  is  10.  Eq  (45)  becomes 


73  * 5 x 10  s i Ji^  to 

~ 5.7  + I.  o 

4.7 


(46) 


Figure  11  shows  that  A/B  exceeds  this  threshold  everywhere 
except  very  close  to  the  two  null  points.  If  either  Tj  or 
the  ratio  is  increased,  the  area  of  acceptable  operation 
will  decrease.  This  example  illustrates  the  procedure  that 
can  be  followed  by  the  engineer  in  evaluating  his  system's 
operation. 

In  the  course  of  studying  the  properties  of  this  cr ray , 
a oomputer  program  was  written  to  implement  Eq  (14).  The 
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appendix  discusses  the  details  of  this  program.  The  only 
important  details  needed  here  are  that  the  parameters  of 
this  example  were  included  in  the  program,  plots  for  various 
signal  locations  were  requested,  and  the  signal  used  was  the 
triangular  signal  of  Figure  12. 

As  is  easily  shown,  xhe  first  derivative  of  a trian- 
gular signal  is  a square  wave.  The  basic  form  of  the  sig- 
nal and  its  derivative  were  illustrated  in  Figure  12.  This 
signal  was  originally  selected  because  it  represents  the 
correlation  function  of  a PN  code.  Notice  that  the  time 
derivative  has  an  amplitude  that  is  5 x 10^  times  larger 
than  the  signal.  This  will  cause  the  time  derivative  term 
to  have  more  effect  on  the  overall  output  than  predicted  by 
the  constants  A and  B alone.  This  should  lead  the  reader 
to  the  realization  that  signals  with  triangular  correlation 
functions  are  not  good  signals  to  use  with  arrays  if  the 
effect  of  derivatives  is  undesirable. 

As  will  be  shown  momentarily,  the  plots  from  this  com- 
puter program  illustrate  quite  nicely  the  results  discussed 
to  this  point.  By  substituting  the  parameters  of  this  ex- 
ample into  the  program  we  are  in  essence  programming  the 
time  domain  version  of  Eq  (37)  which  is  basically  Eq  (39). 

The  computer  program  was  run  for  seven  different  signal 
locationsi  n° , 1°  , 10*  , 40° , 60° , 70.53'’.  and  90# , These 
locations  include  the  two  null  positions,  0°  and  60*,  and 
the  main  beam  axis,  90*.  The  output  waveform  is  complex. 

It  can  be  manipulated  in  several  different  ways  to  obtain 
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the  equivalent  output  of  real  systems. 

The  first  operation  performed  is  to  plot  the  real  part 
of  the  complex  waveform.  This  is  equivalent  to  looking  at 
a coherent  detector  that  is  locked  onto  zero  relative  phase. 
As  can  he  seen  from  the  phase  shift  terms  of  Bq  (39) » the 
phase  of  the  signal  varies  from  2?0°  to  0°  (passing  through 
180®)  as  the  signal  location  moves  from  0®  to  90®  while  the 
phase  of  the  derivative  term  moves  from  180°  to  -90*  in  the 
same  interval.  Thus,  at  0°  and  at  cos“*(l/3)  (approximately 
70.53°)  it  is  seen  that  the  time  derivative's  phase  places 
it  on  the  real  axis  and  the  signal's  phase  places  it  on  tno 
imaginary  axis.  Thus,  this  coherent  receiver  will  not  see 
the  signal  at  these  locations.  At  other  signal  locations, 
the  receiver  will  see  components  of  both  signals  on  the  real 
axis.  The  size  of  the  signal  will  he  large  compared  to  the 
derivative  in  most  casts  and  the  derivative  won't  be  ob- 
servable. Near  60°,  though,  the  signal's  magnitude  drops 
and  the  derivative  component  becomes  observable.  Figures 
13  through  19  illustrate  what  the  real  part  of  the  output 
looks  like.  As  can  be  seen,  the  output  is  mostly  the 
derivative  near  the  two  nulls , 0*  and  60" , and  it  is  mostly 
the  derivative  at  70.53°*  At  10*  we  observe  the  effect  of 
both  signals  on  the  output.  If  we  had  selected  a signal 
closer  to  the  null  at  0® , the  derivative  term  would  be  more 
evident  in  the  output. 

The  reader  should  note  that  the  figures  that  contain 
the  square  wave  signal  have  low  amplitude  oscillations  near 
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the  corners  of  the  square  wave.  These  oscillations  are 
caused  by  the  Gibbs  phenomenon  (Ref  2il05-10?).  In  essence, 
these  oscillations  are  created  any  time  a signal  with  dis- 
continuities is  represented  by  a finite  length  Fourier 
series. 

The  next  operation  performed  by  the  program  is  to  plot 
the  magnitude  of  the  complex  waveform.  This  is  equivalent 
to  looking  at  the  output  of  an  envelope  detector.  The  re- 
sults are  shown  in  Figures  20  through  26.  Comparing  these 
plots  with  Figures  13  through  19  shows  that  the  envelops 
detector  yields  the  desired  signal  as  an  output  more  often 
than  the  coherent  receiver  usedi  e.g.,  compare  Figure  18 
with  Figure  25.  Figure  21  is  the  most  interesting  of  this 
group.  It  illustrates  the  appearance  of  the  output  when 
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the  signal  term  and  the  derivative  term  are  both  about  the 
same  magnitude. 

The  third  operation  performed  by  the  program  is  to  sep- 
arate the  signal  and  its  derivative  into  separate  quadratures. 
This  is  possible  because  the  signal  and  the  derivative  are 
separated  by  90*  of  phase  shift.  Recall  that  this  is  not 
necessarily  true  for  all  arrays.  It  holds  for  all  linear, 
equally  spaced  arrays  and  this  example  is  a member  of  that 
set,  but  it  does  not  hold  for  all  array  geometries.  Thus, 
we  can  separate  the  signal  and  its  derivative  here,  but  they 
cannot  be  separated  into  separate  quadratures  in  any  system. 
The  technique  used  in  the  program  is  to  multiply  the  complex 
output  by  the  complex  conjugate  of  the  phase  of  the  signal 
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term  in  Eq  (39)*  This  multiplier  is  d * 
result  of  this  operation  is  to  leave  the  signal  as  the  real 
part  of  the  complex  waveform  and  the  derivative  term  as  the 
imaginary  part.  The  program  then  plots  these  two  waveforms 
separately.  The  results  are  shown  in  Figures  27  through  40. 

Figures  27  and  35  show  that  the  signal  term  vanishes  at 

the  two  nulls  and  Figure  40  shows  that  the  derivative  term 

vanishes  at  90°  , the  main  beam  axis.  These  figures  quite 

nicely  illustrate  the  effect  of  the  constants  A and  B from 

Table  I on  the  magnitude  of  the  output  terms.  For  example. 

the  signal  triangle  peaks  at  2 x 10**^  (Fig.  12a).  The  para- 

—4  * 

meter  A has  the  value  of  -9*5  x 10  at  1 (Table  I).  Thus, 
the  output  of  the  signal  quadrature  should  be 

2 x 1(T6  x -9.5  x 10~4  = -1.9  x 10"9  (47) 

As  can  be  seen,  this  is  the  peak  value  of  the  triangular 
waveform  in  Figure  29.  This  result  can  be  confirmed  in  all 
the  other  figures  shown  here. 

The  inverted  waveform  of  Figure  29  brings  up  an  inter- 
esting point  about  threshold  detector  systems.  If  this 
array  was  used  with  a system  that  transmitted  positive  going 
pulses  to  represent  a +1  and  negative  going  pulses  to  repre- 
sent 0,  then  the  detector  would  have  interpreted  the  pulse 
at  1®  incorrectly.  Figure  37  shows  that  the  same  pulse 
would  have  been  interpreted  correctly  if  the  signal  were  at 
70* . With  a little  thought  the  reader  should  realize  that 
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this  effect  would  occur  with  any  antipodal  signal  set. 

Thus,  a threshold  detector  of  this  type  cannot  operate 
properly  with  an  array  unless  it  has  some  additional  me- 
chanism  to  determine  the  expected  Bign  of  the  signal. 

The  last  operation  performed  by  the  program  is  to 
calculate  the  phase  of  the  waveform  and  plot  the  results. 
Before  looking  at  the  plots  it  is  helpful  to  obtain  a feel 
for  the  type  of  result  to  expect.  An  example  will  bring 
out  all  the  relevant  ideas. 

To  simplify  the  discussion,  assume  that  the  term  A 
and  B from  Table  I are  such  that  the  output  signal  (the 
triangular  waveform)  and  the  output  derivative  (the  square 
wave)  have  the  same  peak  gain.  Since  these  two  outputs  are 
in  quadrature  the  output  can  be  represented  in  a three  di- 
mensional representation  of  the  form  of  Figure  4l.  The 
triangular  waveform  is  in  the  0°  - t plane.  The  square 
wave  is  in  the  90°  - t plane.  The  system  output  is  then 
the  vector  sum  of  these  two  signals. 

From  this  drawing  it  can  be  seen  that  the  system  out- 
put phase  will  start  at  90°  at  t = -T  . As  t goes  from 
t = -T  toward  t = 0 , the  phase  will  vary  from  90*  to 

45®.  As  t passes  through  t = 0 , the  phase  will  instan- 

taneously switch  from  45°  to  -45°.  Then  as  t proceeds 
toward  t = T , the  phase  will  transition  from  -45*  to 
-90°. 

Although  a real  system  cannot  instantaneously  change 
phase,  this  figure  indicates  that  the  phase  makes  a 
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significant  and  very  rapid  phase  change  over  the  duration  of 
the  signal  when  the  signal  and  its  derivative  are  both 
equally  present,  A phase  tracking  system  would  be  hard 
pressed  to  track  this  degree  of  phase  shifting. 

Now  the  effect  of  changing  the  relative  sizes  of  the 
signals  in  Figure  4l  becomes  more  apparent.  As  the  signal 
term  is  increased,  the  overall  phase  moves  closer  to  the 
signal  phase.  Figures  42  through  48  show  the  results  from 
the  computer  program. 

The  figures  need  some  clarification  before  they  become 
understandable.  Figure  42  is  the  plot  for  0°.  At  this 
point,  the  signal  is  zero  and  the  phase  plot  should  just 
be  the  phase  of  the  derivative  term  which  is  either  90°  or 
-90°  depending  on  its  sign.  However,  the  calculations  made 
by  the  computer  do  not  yield  exactly  zero  results  for  the 
signal  quadrature.  The  calculations  yield  a signal  term 
that  oscillates  above  and  below  zero.  The  computer,  using 
an  algorithm  for  tan”*(x)  that  ranges  between  90°  and  -90°, 
calculates  the  phase  to  be  -89.9°  when  in  reality  it  would 
have  been  90.1° . Thus,  when  the  system  phase  oscillates 
about  90*  or  -90° , the  computer  calculates  the  phase  os- 
cillations to  be  wildly  fluctuating  between  90*  and  -90c . 
With  additional  programming,  this  effect  could  be  elimi- 
nated. However,  it  only  occurs  when  the  phase  is  calcu- 
lated at  a null.  Thus,  it  is  left  to  the  reader  to  under- 
stand the  source  of  these  graphic  fluctuations  at  the  two 
null  points. 
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Figure  42.  Output  Signal  Phase  of  Linear  Array. 
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Figure  45.  Output  Signal  Phase  of  Linear  Array 
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Figure  48.  Output  Signal  Phase  of  Linear  Array. 


Figure  4-3  shows  that  the  phase  changes  dramatically 
in  the  vicinity  of  nulls.  Figures  44,  45,  47,  and  48  show 
that  the  phase  remains  close  to  zero  when  the  signal  term 
is  dominant.  Thus,  in  mopt  regions  a phase  tracking  sys- 
tem can  operate  without  d*-ficulty.  However,  near  the  nulls 
typical  phase  tracking  systems  will  probably  lose  track  due 
to  the  fast  changes  in  the  output  phase . 

The  computer  program  illustrated  the  effects  of  the 
array  on  incoming  signals  by  using  a triangular  signal. 

This  signal  also  yields  insights  into  the  signal/time  deriv- 
ative relationship. 

As  was  mentioned  earlier  this  waveform  was  initially 
selected  because  of  a desire  to  study  the  array  effects  on 
the  correlation  properties  of  PN  code  pulses.  Nhat  resulted 
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is  summarized  in  Figure  49.  Si^oly  stated,  the  figure 
illustrates  that  the  derivativt  a pulse  has  more  effect 
on  the  output  of  the  correlator  than  the  pulse  whenever  the 
pulse  width  is  less  than  one  second.  For  a $ MHz  pulse 
rate,  the  effect  of  the  derivative  term  is  5 x 10^  times 
greater.  This  gain  term  will  offset  a lot  of  the  insig- 
nificance of  B in  Table  I. 

The  basic  results  obtained  from  studying  this  four 
element  array  can  be  summarized  as  follows « (1)  the  array 

equation  doe.,  result  in  the  output  of  both  the  signal  and 
its  time  derivative}  (2)  for  this  example,  the  signal  and 
its  derivative  are  phase  shifted  by  90° i (3)  in  most  areas 
the  signal  output  is  dominant j (4)  in  the  neighborhood  of 
nulls  the  phase  of  the  output  signal  varies  dramatically. 
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D.  Sparse  Arrays 


In  the  previous  section  a computer  program  was  used  to 
illustrate  the  effects  of  a linear,  equally  spaced  array  on 
the  incoming  signal.  The  program  just  illustrated  the  re- 
sults that  had  already  been  developed  analytically.  In  this 
section  we  use  the  computer  program  to  obtain  results  for 
two  arrays  that  cannot  be  analyzed  in  closed  form. 

The  first  array  to  be  analyzed  is  a four  element,  linear 
array.  However,  the  spacing  is  unequal.  The  first  three 
elements  are  spaced  one  half  wavelength  apart  and  the  fourth 
element  is  located  five  wavelengths  from  the  third  element. 
This  type  of  element  arrangement  is  called  a sparse  or 
thinned  array.  The  study  of  sparse  arrays  is  useful  for 
several  reasons.  First,  they  are  still  fairly  easy  to  ana- 
lyze because  of  the  rotational  symmetry  of  the  antenna  pat- 
terns i that  is,  the  output  waveform  is  still  independent  of 
^ (see  Fig.  1).  Second,  in  applications  where  the  antenna 
elements  are  mounted  or  an  aircraft  surface,  the  aircraft 
configuration  may  force  a sparse  arrangement  on  the  designer. 
Third,  and  probably  the  most  important,  is  the  fact  that  the 
elements  are  not  really  isotropic  (due  to  their  ivn  design, 
due  to  the  aircraft  structure).  This  may  force  the  designer 
to  arrange  the  elements  for  a larger  field  of  view.  Finally, 
although  it  is  not  proven  here,  a properly  designed  sparse 
array  can  have  some  of  the  properties  of  larger  arrays  with- 
out tb\.  expense  of  the  additional  elements  (Ref  7» 121-126). 

For  the  purpose  of  this  analysis,  the  weighting 
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coefficients  will  be  assumed  to  be  real  and  of  value  one 
(i.e.,  fi.'t  , Hji  * 0 ).  This  has  the  effect  of  steering 

the  main  beam  to  broadside  (0  = 90* ) . The  analysis  focuses 
on  the  first  two  terms  of  Eq  (25).  This  part  of  the  equa- 
tion can  be  written  as 

Z (6,*)*  A (G)m(i*tc)*  (48) 


where 


A*  I 


(49) 


and 


6(8) 


jit 


C os  & 


./Jtvf,  Cdj  - de)oo » ® 


1 


(50) 


where  di~d-e  - - I f , ^ - cf4  * -|f  ? -dc  s ' ? , and  -de  - Vi . 

The  analysis  begins  by  evaluating  A(0)  and  B(0)  in 
terms  of  0.  Because  of  the  symmetries  involved, 

A(0)  = A(-0)  and  B(0)  = B(-0)  . Furthermore, 

A(90‘  - 0)  = A(90*  + 0)  with  similar  results  holding  for 
B(0).  Thus,  all  values  of  A(6)  and  B(0)  can  be  easily  ob- 
tained from  values  of  A(0)  and  B(0)  in  the  first  quadrant. 
Table  IY  lists  values  for  these  two  parameters  in  both  rec- 
tangular and  polar  format  for  the  first  quadrant.  The  last 
column  of  the  table  also  shows  the  ratio  of  |a)/|b|  for  this 
array. 
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We  recall  from  earlier  analyses  that  A(0)  is  the  usual 
array  factor.  It  is  graphed  in  Figure  50  in  a polar  plot. 
This  would  normally  be  called  the  antenna  pattern  for  this 
array.  Figure  51  is  a polar  plot  of  B(0).  We  called  this 
the  time  derivative  radiation  pattern  in  the  last  section. 
The  first  interesting  thing  to  note  about  Figure  50  is  that 
it  is  vastly  different  from  the  four  element  pattern  of 
Figure  9 1 it  has  a quickly  changing  antenna  pattern  with  a 
lot  of  sidelobes  and  dips  in  the  gain.  The  second  thing  to 
note  is  that  there  are  no  deep  nulls  for  this  array. 

Having  no  deep  nulls  results  in  two  consequences. 

First,  the  signal  term  is  never  completely  attenuated.  This 
is  confirmed  by  the  A(0)  term  in  Table  IV.  Second,  since 
A(0)  is  never  nulled,  it  and  B(0)  are  never  comparable  in 
magnitude.  The  ratio  of  their  magnitudes  is  never  less  than 

O Q 

10  . Thus,  the  derivative  term  will  have  to  be  10  times 
larger  than  the  signal  term  in  order  ior  the  derivative  term 
of  the  output  to  be  equal  to  the  signal  term.  If  we  use  the 
signal  of  Section  III.  C.  on  this  array  and  an  envelope  de- 
tector, the  signal  will  always  be  the  dominant  part  of  the 
output . 

Comparing  Figures  SO  and  51  illustrates  that  A(0)  and 
B(0)  are  still  reciprocally  related  as  was  noted  with  the 
equally  spaced  array  of  the  last  section.  As  was  stated 
earlier,  though,  this  relationship  is  simply  stated  as  an 
observation  here  and  no  attempt  has  been  made  to  analyze 
this  apparent  inverse  relationship. 
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Maximum  Value 


Element  Locations  in 
Wavelengths 
0.  0.5.  1.0.  6 


Figure  50*  Four  Element  Sparse  Array,  Plot  of/A(0)( 


Max.  Value  8 1.4  x 10 


Element  Locations  in 
Wavelengths 

0*5 , 1.0,  6.0 


Figure  51*  Four  Element  Sparse  Array,  Plot  of  [i>(0)| 


In  the  last  section,  the  signal  and  its  first  deriva- 
tive were  always  90°  out  of  phase.  By  looking  at  the  phase 
terms  in  Table  IV , we  see  that  A(6)  and  B(0)  are  not  sepa- 
rated by  90° . In  fact,  they  are  not  even  separated  by  a 
constant  phase.  There  is  only  one  way  to  isolate  only  one 
signal  term  with  this  array.  To  see  the  derivative  term, 
the  system  would  have  to  track  the  phase  of  A(0)  and  then 
coherently  detect  the  output  at  a phase  that  is  90*  removed 
from  the  phase  of  A(0).  The  output  would  then  contain  only 
that  portion  of  the  derivative  term  that  projects  onto  this 
phase . 

If  one  is  interested  in  being  able  to  detect  the  signal 
regardless  of  its  location,  this  sparse  array  is  a good 
array  in  that  it  never  nulls  the  signal  term.  However,  the 
usual  purpose  for  using  an  array  is  to  be  able  to  place 
nulls  on  interference  sources.  Since  this  array  has  no 
nulls,  the  interference  source  cannot  be  nulled. 

At  this  point  we  drop  discussion  of  this  particular 
sparse  array  and  focus  attention  on  the  second  sparse  array. 
This  array  is  identical  to  the  first  array  except  that  the 
fourth  element  is  located  4.5  wavelengths  from  the  third 
element  instead  of  5*  Table  V shows  values  for  A(0)  and 
B(6)  for  this  array.  Figure  52  is  the  polar  plot  of  A(0) 
for  this  array  and  Figure  53  is  the  plot  of  B(0).  They  are 
somewhat  different  from  the  plots  obtained  for  the  first 
sparse  array.  A(8)  still  has  a large  number  of  sidelobes 
and  dips,  but  in  thiu  case  there  are  nulls  associated  with 


< 


82 


mu  v.  *4utt  or  * mn  . re*  »«  t»*mr  »«w», 


miS  PA(ffl  IS  BEST  QUAL1W  mC.« 
UtOM  COPT  f vlfflULSH®  TO 


« 

tie 

A 

ACU 

l«10TN6ftf 

• 

•m 

I*«6Xh8*V 

1*1 

A 

OFCACft 

H 

& 

oriAc** 

i»/ti 

t 

-*l4)r-0> 

•ellTMl 

t9fd?*eo 

OO.I 

*141r-07 

•100.0 

,11«F* 

r 

-*164f-65 

• 18INII 

-•94F*- ll 

• 1U7.fl 

18.8 

• 14  1C*  97 

-153.0 

,l?4f* 

s 

•*1437-04 

•4M7-01 

-•Utir-or 

•9HP-0I 

90.0 

•143C-07 

-150.9 

• n?r* 

* 

-*t«SE»Q4 

•*t4?C-J> 

.7UC-C! 

90.0 

•I*:f-07 

•179,4 

.5377* 

e 

•*i4ir«63 

•tnr*r« 

-.t4:t-e> 

90*1 

.14tV-3> 

-179.9 

.0417* 

» 

-•795C-03 

-.1417-07 

•i7r7*oi 

90.1 

•141E-07 

-170.0 

• rr» 

r 

•.045C-CJ 

•7iu*oe 

-.mr-op 

•%19F*C0 

90.1 

•t4ir-Q? 

-170.1 

•1467* 

i 

•104* *00 

• •Hjr.o? 

-.99*»r-l3 

•3A97*0» 

99.7 

•tftJP-07 

•17). 7 

t 

•tVU!»0P 

• «»  r*oo 

-.U47*0? 

• .i3r»--8ti 

.31^*10 

90*2 

.1 J4F-07 

-17010 

,!7)*» 

i« 

.4>*f*C0 

-. 1367-07 

••177F-39 

.9f::*oe 

90*1 

• 1 lvr- 57 

-170.5 

it 

••ins-»? 

,9M*  *91 

-.1*47-37 

.93!t*08 

90.1 

.1347-97 

-t?S.4 

.4r»r* 

w 

t^,f*81 

-.i ur-or 

-•ltOh-09 

.97;  :*08 

9U9 

.1  JlC-37 

-179.3 

• 614N 

ti 

-.5n?-37 

•701**08 

-.Uic-or 

•oHlr*8l 

••8i*:*p8 

53.4 

•177C-9? 

-174.1 

.61ft?* 

it 

-.??7F-C S 

.8957*83 

-•U?C-0P 

-•*2 9F-39 

•8ll7«C0 

90.1 

•177C-07 

-179.0 

.7NE* 

it 

•tUSf*#? 

• 1 0 1 F * 0 1 

-tU8£«flF 

-.71  tf-81 

•tftlCtOl 

90.8 

• U7?-07 

•170.0 

,4?cr* 

it 

*.t»t?-oi 

• tiu*ct 

-.ttOC-OP 

••'itr-oo 

• 1 1 ? 7 • 0 l 

90.0 

• no*:- J7 

-170.0 

• 1 0 3 r ♦ 

ir 

•.146E-01 

•t?*r*ot 

-•10’foOP 

••COftp-09 

90.7 

.1077-07 

-1 7 4 .0 

.l?3r* 

it 

-.t>C*>01 

.Il?'*81 

••l83»-0f 

•7Jf7»01 

90,7 

• 9 10C-41 

•IM.9 

• 1-  7 1 * 

it 

-.nst-n 

.t<.4?*oi 

-••u?-e> 

••183F-09 

iur*ei 

90.7 

•tJo'-oe 

•179.0 

ft 

-.7£1C-C1 

.lMNflt 

-•>?in-oo 

-•t10r-09 

.18 Jr*0l 

90.0 

•779E-03 

-l?5.7 

.719?* 

it 

,lf9r»|ti 

•.At4t*03 

*.»!?*•» 

•nu*ot 

90,7 

.6O9r-30 

•179.7 

,?74?* 

ft 

-,M6» -Cl 

.iur*oi 

- •9»*0V  * 08 

• AS>r-n 

.1777*81 

90.7 

•40  3r-01 

179.? 

• J’lH* 

ii 

•.195?-3l 

• mNii 

-.n?r-oi 

• Uir-09 

•10*7*01 

90*9 

• 34?r.-00 

177.0 

.5177* 

ft 

••m?*crt 

• lUr*  31 

• UC«--09 

• 1 IlC*  0 t 

95.9 

•707S-99 

170.0 

.943?* 

>t 

tlllFOl 

-.*.71: -31 

.98C*-*0f 

•1107*01 

90.3 

•6777-39 

134.9 

• 70  If  * 

ft 

*t1«r*81 

• 10H**0» 

.9^9r-89 

•1 907*01 

04.0 

•1731-08 

37.5 

.liar* 

tr 

•t*>f*01 

• «9?«--09 

•1077*01 

04.9 

.767E-40 

10.4 

. 69*6* 

ft 

.W-81 

.188f*81 

• Ult-OI 

•105F-8S 

•i8or.*fi 

08.4 

.413F-90 

14.0 

.4307* 

ft 

,5S5*,-9i 

.1T»T*<*1 

«luf*0l 

•1767*81 

04.0 

.3511*81 

13,7 

.360?* 

•It 

tt$^»8! 

•Atir-OA 

*imf-C0 

•i57,:*flt 

£9.0 

•1771-0. 8 

13.3 

**i?i * 

ll 

,»!J*88 

• 1 4*  7 ♦ f f 

•m*-oa 

•t Alc-08 

•1917*61 

09.0 

•747C-1J 

11.8 

.tf CE* 

it 

.uir  *io 

.IWffOl 

•a&ic-oo 

•t78F-08 

•t*l?*6l 

07a 

•0707-00 

11.4 

• t*»  v C* 

ll 

.99»r*00 

• 9.NS-0I 

•13IF-08 

•to?c*ot 

70.3 

. 9447-00 

11.1 

•loot* 

it 

.t$K*88 

.?•;?«■*  no 

• 9b7F-0  4 

•I47F-CA 

•700r*0O 

71.0 

.99-7- <14 

10.7 

.4117* 

it 

.loot  »0l 

.90«7*08 

•971C-00 

•tffr-o* 

• RU7*oo 

99.0 

asss-oi 

10.1 

,5c  6 E * 

it 

• 346«**oo 

.V,*CI 

.1A5C*0t 

•I58r-O0 

•417  7*00 

34,0 

•976C-O0 

1.3 

,4-?7f* 

if 

• UP?  *10 

«.14<r-Ci 

t«i?E-CI 

,t1D**3S 

•1*07*00 

T0.1 

4.t 

•4t«e* 

It 

• A?0f*.lC 

• *?0.  C * 00 

• 8*t?C-00 

•916F-09 

.5107*00 

-39,0 

•044c-00 

0.7 

,004?* 

n 

,440r *88 

-,H’f*3C 

• 7%S6  »Od 

• 9'»?r-09 

-53.7 

•746r-03 

3.7 

,93?7* 

to 

<445?  **0 

-.7517*00 

.A7*t*:-0ft 

•8707*00 

-59.9 

.o:r.fri0 

•4«6rVo8 

-1.5 

• 145C* 

ti 

• 4?<IC*5C 

*.i«r»c8 

•904F-00 

•ti:if*M 

•tan*01 

-99.4 

•9.0 

.284?* 

tt 

,Mr  a 

•.t0«r*11 

•391C-0I 

•ttJl'-08 

ar«r*3i 

-70.3 

•376C-03 

-?9 .0 

,1C4^* 

ti 

-.trr*ot 

• mc*03 

•,?Uf*C4 

.1717*01 

-74.0 

.3)11-10 

•9  0.9 

,4«IF* 

tt 

tnitoo 

-.i?tr*ot 

• AsK-Oi 

-.1107-08 

•1?J7*01 

-79.0 

•3137-34 

-M.9 

,394f* 

tt 

-.t?0r*01 

-.9lte-09 

-.T78r-ca 

•iiir*oi 

•09.0 

•4097-30 

•133.6 

.315F* 

tt 

•tJUe*ll 

•.U*f*M 

-aicc-oo 

•1107*01 

•98*4 

• 40  7f -08 

-115.7 

,?U'* 

ti 

-.♦rt*:  *oo 

••1047*01 

-..U»*r-PO 

-•••fvr-oa 

•1077*01 

•141.0 

• 544f - a 0 

-1.U.7 

. m »■  * 

tt 

-,461»*,3,i 

-.497r*e8 

-.17?f-00 

-.99»*r-e* 

.101* *01 

-117.7 

.6157-00 

-176.9 

■ it 

tt 

-•AN*:  *0  3 

•t»«J*81 

-•^0^£-04 

u*-ti 

•001 ►♦00 

-1  1?.9 

•0167-00 

•UO.O 

• :v>*n 

tt 

• a»4f  *oo 

•.541«-*00 

-.9I0K-J0 

-.ftT*‘,-58 

.1047*01 

•t  4ft  »6 

.6U**-30 

-H?. 5 

,i4t;r* 

ti 

-.U3?«n 

-.iMr*oe 

-.39JC-03’ 

•1147*01 

•19?. 0 

• 95.»t:*oi 

-i:<».7 

,?‘*to7* 

ti 

-.IMF  *91 

-.17*7*00 

-.393C-09 

-.?98r.C8 

•17?. 1 

.•9-7- !li 

-iJi.i 

*2.»ir* 

ti 

•,npf«ci 

-,t7*r-oi 

-•375F-04 

-•OOr-38 

.1  J7f.*0l 

•179.3 

-1-*6.1 

,417V* 

tt 

%i^r*ai 

• 10«  F « 00 

-•1O4C-04 

-•9?lF-09 

175.0 

•19«F -30 

-167.7 

,7'3»* 

tt 

•-.t****ai 

.t^r*co 

-.117E-04 

.1^97-08 

177.0 

.1517-00 

114,7 

.911*-* 

tt 

-.utr*:i 

• ?.*tr*fo 

. J6*tT-P9 

•13*7*oi 

170.4 

•7657-30 

5?.1 

tl'lC* 

tr 

-•me  *o: 

• 974C-0-I 

air^-oi 

•1117*01 

103.1 

• 4 19*>19 

9U 4 

.??  f«'* 

tt 

-*3?ir*co 

•167? *10 

•lc^r-co 

, .59 lr- 00 

.8307*00 

164.9 

.55)7-09 

>9.3 

,tv:r* 

it 

-*446?*O0 

•9t*r-ct 

•1?*E“C0 

.OAff-OS 

•9Ut*C3 

160*4 

• 661 V -0  0 

7$,? 

tin** 

tt 

•isen-is 

-.??=r-ift 

,ujr-co 

. 7 1 9C- 00 

-:t.J 

•7771-O0 

70.7 

,V*5f- 

ti 

,495?  *90 

-.if-'-oi 

•l?r£-9A 

•7i8?-oa 

•*i*-:*03 

-10.9 

.749E-30 

00.1 

tft>?r* 

t! 

• lOtCOl 

-.UCFftCO 

•oeir-oi 

•7197-08 

•191f*fft 

-9.9 

• ?n?-08 

1>.9 

,143*-* 

tl 

• ir>ie  *o  i 

-,?t6!:*ca 

• N5?-09 

• 49  t F-00 

.1307*01 

-a.i 

.6471-08 

49.9 

, 7 10r* 

tt 

• 196C  *C  1 

-.u*f*co 

-•T94C-0* 

.1107*01 

-A. 3 

•C75C-00 

94.7 

.17tf* 

tt 

•?i??*oi 

-*16V*C0 

-.9>5t.-C9 

,1fr.r-0J 

.*117*01 

-4.1 

•381C-O0 

114.5 

.otir* 

tt 

• 750CO1 

-.6Ur-01 

-.itfF-oa 

.Hcr-98 

*3307*01 

-1.4 

.745r-90 

179.3 

*t:?7* 

tf 

t*0kt*n 

*1V?-Cl 

-•70IC-39 

•i^Mr.jo 

t*S47*oi 

7.1 

• 7U7-08 

-:m.i 

tl 

• ?55f  *C  l 

,-®*r*ee 

-.«8t-88 

••1977-98 

.•*377*01 

9.7 

.3011-09 

-119.7 

• 0f-4  7 * 

tt 

,?VE*0t 

•5r*r*co 

-.IMC-0S 

-•17  77-00 

•7117,01 

13.0 

.4341-00 

•\?l *0 

,u  )?•:* 

ft 

oHttOl 

• U-NPtt 

-.ntr-oo 

-•rtl7»-03 

.1107*01 

77.6 

•£34r-30 

-UN? 

,l.Mr* 

fl 

•f t9f«3t 

.<831r*09 

-.lU8p-00 

••670F-08 

•13.11*01 

30.9 

*9471-08 

- 105.1 

*?30C* 

ft 

•%4?E*00 

•tt*7*ci 

-.107C-38 

-•979r-O0 

•1?0C»  01 

60.4 

.0871-00 

-99.1 

•ire:* 

fl 

•tW«8fl 

.171r*0l 

-.354E-89 

-•9797-90 

•1777.pi 

197.7 

•675C-O0 

- >3.9 

*!*>*'* 

ft 

••i?ir*oi 

• l?1r  *«** 

-•9?9f-04 

•17:7*01 

1?5.7 

.6767*04 

•06.4 

,7'6?* 

ft 

-*r 3i;*et 

• liP  f *0t 

ein?-2l 

-•5107-00 

•71 J f*o  t 

191.7 

*5.471-60 

•70.7 

.477** 

ft 

-•?:  jc*oi 

.em*oo 

*l<»9?-08 

-.907F-C8 

•%07 7*01 

167.0 

*t  16!>30 

•07.? 

*6? >7* 

ff 

«,i?5?*oi 

•5M**oo 

•7l3r-08 

-.?9?N08 

*i>t<:*oi 

ICS. 9 

• J3Or-P0 

•50.1 

.1007* 

ft 

•tSf6!»81 

• I9*r  * 90 

•?30?-08 

-• t9?r-C0 

•1577*01 

176.0 

.7511-10 

-?3.9 

,l4?f* 

ft 

*,iMr*oi 

-•5J*?*00 

•??7f-08 

• 9l?‘-09 

*^!rr*oi 

-176.7 

•7311-00 

tO. 3 

• 1571-* 

II 

-,?4Cf *01 

-.Mfr*oo 

*?00f -08 

• t 0 7r- 00 

.1477*01 

-If  9.0 

.7611-90 

15.3 

*119*.* 

tl 

•.?0?£*J  1 

••lUr*91 

• U.or-38 

.^9«r-ce 

• !t>»'*01 

-159.1 

..i3or*o0 

56.4 

.1341* 

It 

•t»ft*9l 

• . t4J  F *0t 

• ltlF-09 

•1017-08 

tfM*  *61 

-144.1 

.1771*90 

64.6 

*t?3C* 

11 

-•lU?*0I 

•tl»fF*M 

• A8*C«0* 

•lirr-oi 

aur»pi 

-174.9 

.3191*08 

>7.0 

•64 1r* 

tt 

-*j7rr*oo 

-.t|7r*81 

.npt-09 

•J89F-80 

•10t?*0t 

-99.3 

.7907-00 

40.1 

.ft6*F» 

tt 

tfircoe 

• .ti<T*oi 

-.If  ftt-39 

• 191 *-C8 

•?04f*0 1 

-67.4 

.7417-30 

94.1 

,1*71  * 

tt 

• ill?  *3  1 

-.tr*c*Pi 

• **tC0?-e<i 

• ir*r-08 

auc*ei 

-4  3.8 

•10?t*30 

t?J  .9 

•i?ir* 

tf 

.»?  )E*.n 

-.14)7*01 

-.9Mtr-c9 

..lttlr-08 

•1071*01 

-74.1 

•1711-04 

113.7 

.ntf* 

tt 

,14??  *01 

•.t**7*Ol 

-aipc-oa 

tirr-;9 

.118.7*01 

-16*9 

.6747-09 

1M.1 

It 

*1357*01 

•tNlf  *00 

39 

.1187-09 

as4«:*f  1 

-8.0 

.7761-0) 

150.1 

• 1*  lr* 

tt 

•931? *01 

*1J»f-i3 

•1357-JI 

•400  E*  0 1 

*0 

.3447. zs 

•a 

•161?* 

17 


'1# 

:a 

H 

1i 

19 

w 

U 

19 

19 

10 

*9 
! 0 

10 

r*» 

5* 

n 

*•> 

19 

J9 

09 

01 

n 

31 

31 

14 

*59 

J1! 

•ft 

09 

‘9 

!9 

:•» 

» 

•“> 

:* 

01 

o* 

09 

;f 

11 

>9 

5* 

0? 


83 


< 


Maximum  Value  = 4,0 


Figure  52.  Four  Element  Sparse  Array,  Plot  of  | A(G)| 


Max.  Value  = 1.4  x 10“^ 


Element  Spacing  in 
Wavelengths 


Figure  53*  Four  Element  Sparse  Array,  Plot  of  |BV6)|. 


it.  Thus,  in  the  neighborhood  of  the  nulls  at  0*  and  60* 
the  derivative  term  should  become  the  dominant  output  of  the 
array.  The  computer  program  was  setup  for  this  array  and 
the  output  at  0*  is  shown  in  Figure  5^*  Ab  is  evident,  the 
derivative  term  comes  through. 

In  the  last  section  we  used  the  ratio  |a(0)/B(0)|  to 
obtain  a measure  of  the  relative  importance  of  the  signal  to 
the  derivative.  The  plot  of  this  term  was  used  in  conjunc- 
tion with  the  ratio  to  determine  when  the  output 

would  have  a significant  time  derivative  component.  Plots 
of  the  ratio  |a(0)/b(0)|  for  these  two  sparse  arrays  are 
included  ae  Figures  55  and  56.  Figure  55  shows  that  m 

O 

would  have  to  be  about  10  in  order  for  the  derivative  to  be 
significant.  Figure  56  shows  that  this  array  has  two  null 
regions  that  would  make  the  derivative  significant.  The 
ratio  ) would  determine  how  far  from  these  two  nulls 

that  the  derivative  remains  significant. 

By  comparing  the  three  plots  of  B(0)  (Figures  10,  51, 
and  53)  we  can  observe  that  B(0)  is  actually  a fairly  con- 
stant function.  Except  for  the  region  near  the  main  beam 
(0  = 90*)  its  values  do  not  fluctuate  that  much,  A(6),  on 
the  other  hand,  is  observed  to  change  amplitude  quite  a bit 
depending  on  the  value  ol  0.  It  is  a possibility  that  this 
constancy  of  B(9)  can  be  of  value  in  the  design  of  systems 
that  rely  on  the  derivative  term.  For  example,  a jammer 
that  is  designed  to  jam  using  the  signal  derivative  can  be 
assured  of  passage  at  a constant  level  except  along  the 
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main  beam  axis. 

By  observing  the  nulls  in  Figures  9 and  52  and  then 
comparing  the  values  in  Tables  1 and  V,  it  becomes  obvious 
that  the  nulls  in  the  sparse  array  are  narrower  than  in  the 
equally  spaced  array.  This  has  the  advantage  of  cutting 
down  on  the  region  where  the  derivative  is  an  effective 
output  signal. 

These  two  examples  pluB  the  linear  array  of  the  last 
section  indicate  how  much  the  basic  array  pattern  can  be 
changed  by  element  spacing.  There  is  another  aspect  of 
arrays  that  these  examples  can  illustrate.  Phase  steering 
of  arrays  (as  discussed  in  the  next  section)  can  have  the 
same  effect  as  moving  the  elements.  By  electronically 
( delaying  the  signal  on  one  of  the  elements,  it  can  effec- 

tively be  moved  away  from  the  other  elements.  However,  the 
more  practical  use  for  phase  steered  arrays  will  be  dis- 
cussed in  the  next  section. 

The  basic  results  to  be  obtained  from  this  section  are 

i 

summarized  as  follows i sparse  arrays  have  nulls  and  the 
output  in  the  neighborhood  of  the  nulls  will  be  dominated  by 
the  derivative  term.  The  signal  and  derivative  terms  are  at 
different  phases.  Unlike  the  linear  arrays  studied  earlier, 

i 

this  phase  difference  varies  depending  on  the  signal  loca- 
tion. 

E.  Phase-Steered  Array 

There  is  one  last  example  to  be  investigated  before 
leaving  this  chapter  on  static  arrays.  The  main  beam  of 
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arrays  can  be  moved  electronically  through  a technique 
called  phase  steering.  This  i3  done  by  effectively  de- 
laying the  signal  coming  out  of  each  element  of  the  array 
by  a specific  predetermined  amount.  The  array  model  of 
Chapter  II  does  not  explicitly  illustrate  this  particular 
technique.  However,  phase  steering  can  be  viewed  as  a 
linear  phase  shift  in  the  frequency  domain.  This  means 
that  a phase  steered  antenna  can  be  modeled  by  replacing 
in  Eq  (14)  by  the  expression  - J *( f ♦ )4ol  where 
j - r-i  r°s  and  is  the  desired  location  of  the 

main  beam  axis.  This  expression  allows  for  both  a linear 
phase  shift  and  a constant  phase  shift.  With  this  change, 
Eq  (14)  becomes 


o(e,  <p,  P)  * tW)£  e 


(51) 


By  factoring  out  the  phase  center  ter  os  and  then  expanding 
the  summation  U3ing  the  Taylor  series  we  obtain 


-4«  • *t)  r 


C j^rVWJV..  * ...I 


(52) 


Inverse  transforming  yields 
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(53) 


As  can  be  seen  this  is  precisely  Eq  (25)  with  A; 
replaced  by  Ax  - Ati  . Thus,  all  the  results  obtained 
from  Eq  (25)  are  still  valid  in  this  case.  The  signal  and 
its  derivatives  are  still  present.  The  derivatives  will 
still  tend  to  be  insignificant  compared  to  the  signal  due 
to  the  oresence  of  the  factors  ( t - jL  . - jL.  )**  in  the  sum- 
mations.  The  summations  are  still  just  complex  numbers 
that  tend  to  phase  shift  the  signal  and  its  ierivatives  by 
different  amounts.  However,  this  phase  shift  has  become  a 
function  of  and  can  now  introduce  a phase-shift  that 

varies  with  the  steering  term. 

In  the  linear  array  the  signal  and  its  derivative  were 
separated  by  $>0# . To  observe  the  effects  of  the  phase 
steering  term  on  linear  arrays  substitute  the  simple  four 
element  array  parameters  of  Section  III.  C.  into  Eq  (51)* 
It  becomes 


6(e,  f)  * Mf 


(5*0 
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where 


A-  I 


Jc 


Cos  6 


&a  B dss>f<i  an^le. 

Upon  performing  the  summation 

Z (0,  *)  - M(f)  e*p[-j  ^r(c«©-c« e )]  e^[^(co5  0 - cos©.)] 

5jn[jtr  (I^J(cm9  - <U>&  ©,  )] 
a>n[s(l*ft)(cos0  -Cos  9)] 

where  the  two  exponentials  are  the  phase-center  terms. 
By  following  the  example  of  Eq  (4i)  the  trigonometric 
term  is  expanded  as 


(56) 


cos  £.)] 

an 

[tfj 

Cos  ©.)] 

sin 
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where 


• t ( ows  6 - co5  ) 


(57) 


Thus, 


and 


A(e,e.) 


(58) 


Jsm(j  Qco>  («?£»)  ~ 1 

“^TiC) 


(59) 


At  this  point  it  should  be  obvious  that  phase  steering 
is  not  going  to  cause  any  new  effects.  The  signal  and  its 
derivative  will  still  be  90*  out  of  phase.  The  signal  is 
still  dominant  except  in  the  neighborhood  of  the  nulls. 

The  steering  will  just  change  the  location  of  the  main  beam, 
the  nulls,  and  the  sidelobes.  In  fact  it  will  change  all 
the  quantitative  results  somewhat  but  not  any  of  the  qual- 
itative results.  Thus,  it  appears  that  phase  steering  will 
not  cause  any  new  effects  to  the  incoming  signal. 

The  one  area  of  quantitative  interest  in  this  example 
would  be  to  investigate  how  much  the  width  of  the  nulls  are 
affected  by  phase  steering.  In  the  previous  sections  we 
showed  how  the  derivative  term  is  dominant  for  some  small 
angular  distance  around  the  null  points.  It  would  be  of 
interest  to  determine  if  this  angular  distance  varies  much 


due  to  the  steering.  The  author  feels  that  the  effect  is 
small  but  did  not  h«^e  the  opportunity  to  investigate  this 
point  thoroughly. 

F.  Summary  of  Static  Array  B?gu.l.t§ 

This  chapter  has  focused  attention  on  the  effects  of 
antenna  arrays  on  incoming  signals.  Although  we  focused 
our  examples  on  the  effects  to  the  autocorrelation  function 
of  a spread  spectrum  signal,  the  main  analysis  is  valid  for 
any  signal. 

The  truly  unique  result  obtained  from  the  foregoing 
analysis  was  the  fact  that  the  output  of  antenna  arrays  is, 
in  fact,  a combination  of  the  incoming  signal  and  all  its 
time  derivatives.  The  development  that  led  to  this  result 
also  indicated  that  the  signal  will  be  the  only  significant 
output  except  in  the  neighborhood  of  nulls.  In  the  neigh- 
borhood of  the  nulls  the  first  derivative  will  also  become 
significant  due  to  the  attenuation  of  the  signal.  The 
analysis  also  showed  that  the  system  output  phase  changes 
rapidly  near  nulls.  Finally,  it  was  shown  that  the  use  of 
pulses  as  signals  and  correlators  as  detectors  will  accen- 
tuate the  effects  of  the  derivative  and  expand  the  neigh- 
borhood around  the  nulls  where  the  effects  are  observable. 
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The  majority  of  the  effort  in  this  thesis  has  concerned 
itself  with  the  effects  of  fixed  arrays  on  signals.  How- 
ever, most  large  arrays  in  use  today  are  adaptive  in  nature. 
Thus,  we  felt  that  it  was  necessary  to  look  at  aspects  of 
adaptive  arrays  and  analyze  their  effect  on  the  incoming 
signal . 

Adaptive  arrays  use  a variety  of  algorithms  to  control 
the  weighting  coefficients  of  the  antenna  elements.  These 
algorithms  process  the  incoming  signals,  noise,  and  inter- 
ference to  determine  the  values  of  the  coefficients.  If 
these  external  conditions  are  unknown,  then  the  adaptive 
effects  can  be  viewed  as  a form  of  noise.  This  approach 
will  be  taken  in  this  chapter. 

The  array  model  of  Eq  (10)  incorporated  the  adaptive 
nature  of  the  arrays  by  making  the  weighting  coefficients 
functions  of  time.  However,  Chapter  III  assumed  that  these 
coefficients  were  constants.  The  resulting  fixed  array  cal- 
culations yielded  very  important  properties  about  arrays  and 
their  effect  on  wideband  signals.  In  this  chapter,  we  will 
assume  that  the  adaptive  array  is  trying  to  seek  a desired 
static  (steady-state)  condition  and  the  environment  is 
causing  perturbations  around  this  desired  state. 

Eq  (10)  is  used  to  develop  some  results  that  illustrate 
the  effect  that  time  varying,  weighting  coefficients  have  on 
the  signal.  Two  models  for  the  weighting  coefficients  are 
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introduced  and  the  analysis  performed.  It  is  hoped  that 
these  simple  models  will  yield  insight  into  more  complex 
situations. 

To  start  the  analysis,  Eq  (10)  is  repeated  here. 


a («,*>,<) 


v 

E 

«.  i 


A (i)  v 


(10) 


Recall  that  this  is  the  complex  baseband  model  of  the  array 
and  its  output  becomes  the  input  to  the  receiver  processor. 
Next,  recall  that  in  Section  III.  E.,  an  equation  for  phase 
steering  was  developed.  By  combining  those  results  with 
Eq  (10)  we  obtain  a generalised  adaptive  array  equation  that 
includes  phase  steering. 


4. 1 


(bo) 


A.  Sinusoidal  Perturbations 

Assume , 


A4U)e. 


•.«) 


ft  et  ******  * ~ > 


(61) 


where  is  a perturbation  of  the  steady  state 

phase  steering  coefficients.  This  model  can  be  used  for 
several  purposes.  (1)  can  be  considered  a 

noise  jitter  and  the  effect  of  this  noise  on  the  signal  can 
be  calculated  for  a variety  of  noise  levels,  , and  noise 
jitter  rates,  ^ . (2)  Another  approach  would  be  to  assume 
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that  the  are  a modulation  resulting  from  a 

jamming  attempt.  In  either  case  we  are  interested  in  the 
effect  this  sinusoidal  perturbation  has  on  the  incoming 
signal.  By  substituting  these  assumptions  into  Eq  (60), 
the  output  becomes 

S(*f,  4>,i)  s £ A,  • >#<i ) 

*•! 


’VO  1 * m* 

Cm  V 


(62) 


a6  *•« 

The  c?  terms  can  be  expanded  in  a Fourier 

series  (Ref  21ill6)  yielulng 


(63) 


It  -mm 


where  the  are  Bessel  functions  of  the  first  kind  and 

order  1 . Substituting  these  results  into  Eq  (62)  yields 

a(e,<p,  *)  *£  v(i . > -^)e  '-C# 


* ** 


Notice  that  the  i ■ C term  is  the  static  phase  steered 


9? 


(64) 
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array  of  Section  III.  E.,  except  for  the  Jtfg)  amplitude 
weight.  The  .{  , ^ term  ia 

H 

£*(***.■4.)  e>^'( 4 * v •<..)  * J W li.  )J 

*»  i 


(65) 

As  can  be  seen  this  is  the  phase  steered  equation  modified 
by  the  Jj($)  terms  and  then  modulated  by  a carrier  fre- 
quency of 

Thus,  the  noise  jitter  has  caused  two  effects.  The  /$ 
have  changed  the  phase  steered  equation  by  multiplying  each 
output  by  ) . If  the  ,3  are  quite  different,  this 

could  affect  the  basic  array  pattern  immensely.  This  effect 
cannot  be  determined  without  more  information  about  the  -6- 
and  the  array  structure.  If  the  are  all  about  the  same, 
then  they  act  basically  as  an  attenuation  of  the  array  out- 
put. The  term  causes  an  infinite  series  of  modulated 
versions  of  the  signal  to  appear  at  the  output.  Depending 
on  the  modulation  frequency,  , and  the  signal  bandwidth 
these  terms  can  causo  signal  aliasing  effects. 

The  extent  that  these  modulation  terms  interfere  de- 
pends on  both  and  on  determines  how  many  of 

the  Bessel  series  terms  are  located  in  the  bandwidth  region 
of  the  signal.  The  determine  if  the  terms  that  are 

in  the  aliasing  region  are  significant  enough  in  amplitude 
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(compared  to  the  i ■ 0 term)  to  cause  diatortion.  A 
•tudy  cf  the  Bessel  function  will  yie*~  soma  information 
on  this  point. 

Tha  Basaal  functions  can  ba  written  as  powar  sarias 
(Ref  5tl43) . 


.r  i)*"* 


(66) 


Prom  this  it  is  aasy  to  see  that 

J.IO)  * I 

J |(o)  ••  o,  f.,  (6?) 

ThUB,  for  4 ’ 0 * Eq  (64)  reduces  precisely  to  the 

static  phase- steered  array  of  Section  III.  E.  For  sma.ll 
£ (lass  than  .2),  C£(4)».*t  and  Jj (^)s.i  . The 

higher  order  Bessel  terms  are  also  small.  Thus,  for  small 
tha  output  is  mostly  the  phase-steered  output  with  a 
small  amount  of  the  modulated  terms  being  passed.  This 
indicates  that  the  aliasing  effects  should  be  small  and  not 
causa  much  signal  interference.  However,  if  the  become 
larger  than  *2,  the  J {fi ) increase  rapidly.  At  this 
point,  the  effect  of  the  first  modulation  term  will  become 
significant  and  the  value  of  will  become  important. 
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Bq  (64)  is  ths  tins  domain  output  of  tho  array.  If  ws 
psrform  a Fourier  transform  of  this  equation,  the  result  is 


o(«,4>,0«£  if  * A‘}j  -Jfm)  (68) 

i -«•  •'> 


where 

M(0  s-  ^(v(i)  (69) 

In  this  form  the  discussion  above  is  easily  confirmed. 
First,  it  is  simple  to  see  that  the  inner  sum  is  the  static 
phase-steered  array  of  Chapter  III  except  for  the  and 

terms.  It  is  obvious  from  this  form  that  the  3Jj(4«) 
effects  are  equivalent  to  changing  the  weighting  coeffi- 
cients, fiA  and  the  are  equivalent  to  changing  the  «*,  . 

Thus,  if  the  and  take  on  a wide  range  of  values, 

the  array  pattern  will  change  considerably. 

Next,  the  effect  of  the  perturbation  frequency,  ^ , 
is  to  shift  the  inner  sum  up  and  down  in  the  frequency  spec- 
trum. Figure  57  portrays  this  frequency  viewpoint  graph- 
ically. It  shown  that  the  spectrum  at  each  harmonic  is 
different  than  the  actual  signal  spectrum  due  to  the 
and  If  modifications  to  the  static  array  equation. 

The  above  analysis  has  crudely  shown  the  effect  of  a 
sinusoidal  oscillation  of  the  array's  weighting  coefficients 
on  a desired  signal.  Roughly,  it  illustrated  that  the  array 
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Figure  57.  Frequency  Domain  Output  of  Array. 


101 


c»» 


t3EES&S3£afi2^&. 


function  output  is  first  modified  by  the  J and  and 
then  modulated  by  the  harmonics  of  the  sinusoidal  frequency, 
V . It  appears  that  these  effect  are  minor  for  small  per- 
turbation  levels.  It  was  not  shown  what  effects  the  Bessel 
functions^)  and  ^ have  on  the  overall  array  pattern.  If 
the  are  quite  different,  then  the  array  pattern  can  be 
changed  extensively.  This  could  cause  the  same  kind  of 
effects  as  phase  steering  (e.g.,  change  the  pattern  so  that 
the  signal  is  in  a null). 

B.  Gaussian  Noise 

In  this  section  we  assume  a different  model  for  the 
coefficients.  Assume 


* •» 


(70) 


where  mA(*)  are  gaussian  random  processes,  with  means  w,  (*) 
and  covariances  Q,  (7, , ) and  cross  covariances  C Ai,  t tA) 

9 

Eq  (60)  becomes 


*i(0,  <P,*) 


(71) 


where  Jf, (&,  f,  it)  is  defined  as  shown.  Now  we  determine  the 
statistics  of  the  output. 

vie  use  the  accepted  notation  for  expected  values,  £[x]  , 


£[x]  -f  *f(K)c/x 
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to  stand  for 


where  ft*)  is  the  probability  density  function  of  x 
(Ref  3«229).  First,  the  expected  value  of  the  output  of 
Eq  (71)  becomes 


E[8(8.<P,4  *£  *U)  E[er''“)] 


(73) 


The  characteristic  function  of  a random  variable  is 
defined  as  (Ref  3i4l9) 

^(v)  = €['.*''*]  (74) 

Furthermore,  the  characteristic  function  of  a gaussian 
random  variable  it  (Ref  3«420) 

j>x(»)  = exrfjvm  ~ va  x ] (75) 

where  m is  the  mean  of  x and  <r * is  the  variance  of  x.  If 
v * 1 , then  Eq  (74)  reduces  to 

tf>x(l)  " (76) 

This  is  precisely  the  form  of  the  expectation  in  Eq  (73) 
except  for  the  time  dependence.  Therefore,  it  can  be 
evaluated  by  determining  its  characteristic  function  and 
then  letting  v B 1 . This  can  be  accomplished  in  the 

following  fashion.  Let  i t*ik  in  Eq  (73)  • Then 

becomes  a gaussian  random  variable.  As  such,  the 


expectation  can  be  evaluated  using  the  results  of  Eqs  (?4) 
and  (75).  Finally,  by  letting  v * 1 , Eq  (73)  becomes 


E[o(e,<p,  jt)]  ,£  - -<’r‘^‘>  j 


(7?) 


where  **(<**)  s C4(  ^ ^ ) . This  operation  is  valid  for 

any  value  of  ^ , Therefore,  the  result  can  be  generalised 
to 


N 

r[d(0}<P,-<)J  a£  ^(0,  e*f>[jrr>A(t)  - ] 


(78) 


This  shows  that  the  effect  of  gaussian  noise  on  the 
expected  output  is  to  introduce  a phase  shift  and  an  atten- 
uation of  the  output  of  each  element.  This  will  change  the 
antenna  pattern.  The  new  pattern  will  then  determine  the 
effect  on  the  signal.  If  the  nr  are  all  approximately 
equal,  the  phase  shift  introduced  will  have  no  new  effects 
on  the  signal.  Notice  that  the  expected  output  decreases 
in  amplitude  as  the  variance,  cr*  , increases.  In  essence, 
this  says  that  as  the  noise  variance  increases,  the  output 
looks  less  like  the  signal  and  more  like  zero  mean  noise. 
This  should  agree  with  intuitive  notions  about  system 
operations  (i.e.,  the  noisier  the  system,  the  poorer  the 
performance) . 

A final  point  needs  to  be  made  about  Eq  (78).  If  the 
noise  is  zero  mean  and  has  small  variance,  the  expected 
output  of  the  array  is  expected  to  be  the  phase-steered 
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output.  This  is  also  intuitive  in  that  it  simply  says  that 
if  there  is  little  noise,  then  the  system  is  basically  fixed 
and  the  fixed  array  analysis  is  valid. 

The  expected  value  of  the  output  is  important  but  we 
also  need  to  have  a feel  for  how  close  to  the  expected  out- 
put the  actual  output  sample  functions  will  be.  This  leads 
to  a determination  of  the  variance  of  the  output.  The  var- 
iance can  be  calculated  from 

<z\e, <M)‘  E[s(0. 4>>  *) o *(e, 4 - c[ 8& <p, 4 E *[*(<>.  *,*)]  (79) 

where  #■  denotes  the  complex  conjugate.  The  first  term  on 
the  right  side  of  the  eqinality  is  evaluated  as  follows 


If  the  /n.u) .mjU)  are  considered  to  be  jointly  normal  pro-- 
A ' ( 

cesses,  then  the  random  variables  obtained  by  sampling  them 
at  are  jointly  normal.  The  characteristic  function 


of  these  random  variables  becomes  (Ref  3 * 503) 


?*(<)''/}]  (8i) 

By  observing  that  this  result  holds  for  any  jtK  and  by 
setting  v. « i and  v^-i  , the  value  of  the  expectation 
in  Eq  (80)  can  be  determined.  By  doing  this,  Eq  (80)  be- 
comes 


E[t(e, <*>,*)]  £[o*(©  <p,*)] 


M Xi(°,<P,  *)  $*(0,  <P,  *) ex p[j( <nU)  - 

9 

•}'(*)}] 


(82) 


Next,  we  realize  that  the  last  term  of  Eq  (79)  is  just  the 
magnitude  of  Eq  (78).  Thus, 


£ i *i  (6t  <p,  t)  <?,  *)  ' m 3(*^ 


- Si 


+ OTJ(*) 


f-°J— J 


(83) 
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Combining  Eqs  (82)  and  (83),  we  finally  obtain  the  variance 


£ &Ae>0>  *)  */(*,?, 


[e»P  {c.j.<V)j  ''] 


(84) 


The  variance  is  a complicated  function  of  the  signal,  array, 
and  noise  parameters. 

There  are  several  things  that  can  be  noted.  First,  if 
there  is  no  noise,  then  the  (1-^(4,/.)—*  O and  Eq  (84) 
reduces  to  zero  as  expected.  Second,  the  means,  rr^ ( j d),  of 
the  noise  cause  a relative  phase  shift  between  the  cross 
component  signals t i.e.,  the  mean  is  a steering  term  that 
changes  the  variance.  Third,  the  individual  variances  cause 
an  attenuation  of  the  crossterm  products.  Th*s  seems  to  say 
that  the  system  variance  decreases  as  the  individual  vari- 
ances increase.  In  fact,  it  appears  that  the  system  vari- 
ance decreases  to  zero  as  the  individual  term  variances 
approach  infinity.  However,  look  at  the  x - j-  terms.  They 
are  of  the  form 


exploit)]  (ex/>[<r  4(i)]  -i) 


(85) 
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In  this  form  it  is  obvious  that  these  terms  become  larger 
as  the  become  larger.  ThuB,  due  to  these  terms  the 

system  variance  does  not  go  to  zero  as  it  first  appeared. 

As  a final  look  at  the  variance,  assume  the  noise  pro- 
cesses have  equal  means  and  are  uncorrelated.  Thus, 

0 for  i / j . Under  these  conditions 
Eq  (84)  becomes 

<£,  U(e,<*>.  i)|J  (/  -e*p(-<r *<*))} 

4*  I 

4»» 

Thus,  under  these  conditions  we  have  a bound  on  the  system 
variance i a bound  that  is  independent  of  the  noise  levels. 
This  tells  us  that  the  system  variance  level  can  be  com- 
pletely controlled  by  the  designer.  However,  it  is  con- 
trolled by  the  same  parameters  that  determine  the  signal 
level.  Thus,  as  the  variance  is  decreased  so  is  the  signal. 

A measure  of  performance  in  this  case  would  be  to  look 
at  a version  of  the  signal-to-noise  (SNR)  ratio.  Thus,  for 
this  example  of  equal  mean,  independent  gaussian  processes 
the  SNR  ratio  becomes 


snRj  r[g(«.o.<)kf8,(e,«’..o7 

<p,i) 


'-mi] 


t 


t/o  ?,<*)! 4 { 1 

108 


(67) 
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If,  in  addition,  we  assume  equal  variance  terms,  then 


snR*  — — " r it  — 1 — ; 

• •I 


(88) 


The  double  sum  c&i.  be  rewritten  as  the  product  of  the  sums 


SNR  » 


. t *.(«■ P.*)Z  *,'(<>.*■*) 

1 , fl  £ 


>*» 


®*,*(4*1  z 


4»  I 


(89) 


In  this  form  we  notice  that  the  SNR  is  a product  of  a 
noise  controlled  gain  and  an  array  controlled  gain.  The 
noise  factor  shows  that  the  SNR  becomes  very  large  as  the 
noise  variance  approaches  zero. 

Eq  (89)  has  shown  that  the  noise  variance  affects  the 
signal -to-noise  ratio.  Eq  (86)  showed  that  the  noise  var- 
iance is  bounded.  These  two  equations  allow  the  system 
designer  to  pick  the  parameters  to  minimize  the  effects  of 
a jammer  that  somehow  has  managed  to  inject  a noise  signal 
into  the  weighting  coefficient  control  loops.  Prom  the 
jammer' 8 viewpoint,  he  can  use  these  equations  to  determine 
how  much  jamming  effect  he  can  have  if  he  can  inject  noise 
into  the  coefficient  loops. 
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The  analysis  using  sinusoidal  weighting  coefficient 
noist  and  gaussian  noisa  has  revealed  a method  for  degrading 
ths  output  of  ths  array.  This  method  involves  developing 
some  technique  for  affecting  the  control  loops  of  the  array 
so  that  they  are  not  allowed  to  stabilise.  This  implies 
some  very  sophisticated  knowledge  of  the  array  operation. 

The  more  traditional  approach  to  jamming  is  to  simply 
overload  the  system  with  a high  power  interference  source. 
The  standard  measure  of  effectiveness  of  this  method  is  to 
determine  the  signal-to- jammer  power  ratio.  Recall  that 
jf  (6,<P,Jt)  is  the  incoming  signal.  However,  it  can  be 
viewed  as  the  sum  of  the  desired  signal  ar*d  the  jamming 
signal . Thus , 


K(e. <f.  *)  * 


(90) 


where  1^(6,  f,*)  is  ths  signal  component  and  t)  is 

ths  jammer  component.  With  these  definitions  we  can  form 
the  usual  signal- to- jammer  (SJK)  ratio. 


SJ  ft  * 


t *,.(». 4.*) 

ill  * 


Iw'  r 


(91) 


This  ratio  affects  the  overall  system  SNR.  It  shows  very 
dearly  that  one  effective  jamming  technique  is  to  increase 
the  jamming  output  of  the  array.  It  is  not,  however,  as 
simple  as  increasing  the  jammer  power  (though  increased 
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jamming  power  works).  Arrays  are  designod  to  counter  the 
effect  of  the  jammer  by  making  this  part  of  the  output 
small.  The  array  accomplishes  this  by  placing  a null  at 
the  location  (0^,^).  Ir»  Chapter  III  we  showed  that  one 
technique  for  transmitting  through  a null  is  to  transmit 
a signal  with  high  derivative  content.  Thus,  rather  than 
increase  the  jammer  power,  pick  a jammer  signal  with  most 
of  the  power  in  the  derivative. 

The  next  point  to  make  about  this  SJR  is  that  with  the 
use  of  spread  spectrum,  the  jammer  will  also  have  to  counter 
those  anti-jamming  techniques  in  addition  to  getting  through 
the  array. 

In  this  chapter  we  derived  \he  effects  of  two  noise 
models  on  the  output  of  the  array.  Tne  basic  result  was 
that  noise  in  the  array  weighting  coefficients  results  in 
degraded  signal  output.  It  v.as  also  stated  that  a possible 
jamming  technique  is  to  induce  noise  in  the  array  control 
loops.  No  method  was  illustrated  to  show  how  a jammer  might 
manage  to  get  control  of  the  array's  control  loops. 


Ill 
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V Conclusions  &nd  Recommendations 


This  thesis  has  investigated  the  effects  of  antenna 
arrays  on  wideband  signals.  The  effort  began  by  developing 
a complex  baseband  model  for  the  array's  output  that  in- 
cludes the  antenna's  frequency  effects  (Eq  (10)).,  Chapter 
III  converted  this  to  the  equation  for  a static  array 
(Bq  (12)  and  Bq  (14)).  At  this  point  the  array  was  viewed 
as  a filter  and  its  transfer  function  was  obtained  from 
Bq  (14). 

The  transfer  function  was  expanded  in  a Taylor  Beries,. 
This  expansion  was  then  u^ed  to  determine  the  properties  of 
the  resulting  output  waveform.  The  basic  conclusions  drawn 
from  this  investigation  follow. 

(1)  The  output  of  an  array  consists  of  an  infinite  sum  of 
components.  The  components  consist  of  the  signal  and  all 
the  n'th  order  time  derivatives  of  the  signal. 

(2)  Each  of  the  components  is  weighted  by  a complex  number 
that  adjusts  its  phase  and  amplitude  before  the  components 
are  added  together. 

(3)  In  general,  the  coefficient  of  the  signal  term  is  many 
orders  of  magnitude  larger  than  the  first  derivative  term 
coefficient,  the  first  derivative  coefficient  is  many  orders 
of  magnitude  larger  than  the  second  derivative  coefficient, 
etc.  Because  of  this,  the  signal  term  is  usually  dominant. 
The  only  exception  is  in  the  neighborhood  of  nulls,  where 
the  signal  coefficient  vanishes.  Then  the  first  derivative 
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term  becomes  the  main  output  of  the  array. 

(4)  For  wideband  signals  the  derivative  term  tendu  to 
be  orders  of  magnitude  larger  than  the  signal  term.  This 
tends  to  offset  the  small  coefficient  of  the  derivative 
term.  This  means  that  the  derivative  term  becomes  observ- 
able farther  from  the  null  points  than  would  otherwise  be 
expected. 

(5)  The  signal  and  the  derivative  terms  have  different 
phases.  In  the  case  of  linear  arrays  the  signal  and  its 
first  derivative  are  separated  in  phase  by  90° . This  can 
be  used  by  coherent  receivers  to  isolate  either  the  signal 
or  the  derivative.  For  other  arrays,  the  phase  difference 
is  not  fixed  and  in  general  is  not  90° . 

(6)  The  cases  studied  indicated  that  an  envelope  detec- 
tor tends  to  yield  the  signal  as  the  system  output  more 
often  than  cohet<int  receivers.  However,  we  did  not  inves- 
tigate the  effect  an  envelope  detector  would  have  on  the 
PSK  version  of  a spread  spectrum  code. 

Chapter  IV  looked  at  two  representations  of  the  array 
where  the  weighting  coefficients  were  modeled  as  static 
components  plus  an  added  noise  signal.  The  first  repre- 
sentation modeled  the  noise  term  as  a precise  sinusoidal. 
The  conclusions  drawn  from  this  model  follow. 

(7)  The  array  output  is  carrier  modulated  by  all  of  the 
sinusoid  harmonics.  Thus,  the  frequency  spectrum  of  the 
array  output  is  moved  up  and  down  the  frequency  axis  by 
integer  amounts  of  the  sinusoidal  frequency.  If  this 
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f frequency  is  not  high  enough  aliasing  will  occur. 

(8)  Bach  of  tha  modulated  array  outputs  is  modified  by 

\ 

. Bessel  function  coefficients.  These  coefficients  alter  the 

i 

array  output.  Thus,  the  signal  at  one  harmonic  is  not  the 

‘ same  as  the  signal  at  other  harmonics. 

The  second  approach  to  the  study  of  noise  effects 

' assumed  that  the  noise  terms  were  gausaian  random  processes. 

\ 

This  analysis  resulted  in  the  following  conclusions. 

(9)  The  noise  tends  to  drive  the  expected  output  of  the 
array  away  from  the  static  array  output  toward  aero  output. 

(10)  The  variance  of  the  system  output  increases  as  the 
noise  level  increases. 

A crude  look  at  using  noise  to  jam  adaptive  arrays  was 
^ investigated.  One  conclusion  was  drawn. 

(11)  The  jammer  has  two  ways  available  to  change  the  ex- 
pected output  of  the  array.  First,  the  jammer  can  simply 
increase  the  power  of  its  jamming  signal.  This  will  increase 
the  variance  of  the  noise  output  of  the  array.  Second,  the 
jammer  can  transmit  a signal  that  controls  the  variance  of 
the  steering  coefficients. 

The  analysis  performed  in  this  thesis  brought  to  light 
a variety  of  information  about  arrays.  This  information  in- 
dicates areas  where  further  research  might  be  able  to  take 
advantage  of  the  results  obtained  in  this  thesis.  We 
recommend  the  following. 

(1)  The  first  derivative  term  exists  at  the  output  of  the 
array.  Research  should  be  performed  to  see  if  receivers  can 
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be  designed  to  take  advantage  of  tne  presence  of  this  signal. 

(2)  Some  systems  may  want  to  minimise  the  effects  of  the 
derivative  terms.  Research  should  be  directed  toward  find- 
ing signal  sets  that  have  small  derivatives.  More  speci- 
fically, research  should  be  directed  toward  finding  signal 
sets  whose  correlation  functions  have  small  derivatives. 

(3)  For  the  purposes  of  jamming,  signals  wi.h  high  power 
in  the  derivative  terms  would  be  difficult  to  null.  Studies 
in  this  area  could  result  in  valuable  information  about  jam- 
ming of  arrays. 

(4)  One  of  the  original  purposes  of  this  thesis  was  to  in- 
vestigate the  effects  of  the  array  on  spread  spectrum  sig- 
nals. The  details  of  this  analysis  were  never  carried  out 
for  the  specific  spread  spectrum  model.  This  analysis 
should  be  completed. 

(5)  The  effects  of  specific  models  of  adaptive  arrays  need 
to  be  investigated.  Specifically,  transient  effects  should 
be  investigated  when  the  array  adapts  from  one  steady  state 
to  another. 

(6)  It  was  observed  that  there  appears  to  be  an  inverse 
relationship  between  A,  the  multiplier  of  the  signal,  and 
B,  the  multiplier  of  the  derivative.  This  might  imply  that 
the  deeper  the  null,  the  higher  the  derivative.  This  rela- 
tionship should  be  studied  in  detail.  It  might  reveal  in- 
sights into  more  optimum  nulling  techniques. 
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Appendix  t Analysis  of  Computer  Program 


Introduction 

The  purpose  of  this  Appendix  is  to  present  the  details 
of  the  computer  program  used  in  this  thesis.  It  is  divided 
into  four  parts i theoretical  analysis,  computer  program 
description,  discussion  of  data,  and  summary. 

Theoretical  Analysis 

Eq  (14)  is  the  starting  point  with  / and  * 0 . 

Next,  the  array  is  assumed  to  be  linear.  The  lx  become 

~ cos  9 . With  these  substitutions  Eq  (14)  becomes 

4 c 


o(e,  C)  * Hfflle'  c 

I 


(92) 


The  signal  is  assumed  to  be  a pulse  of  unit  height  and  dura- 
tion, T. 

This  waveform  is  passed  through  the  array  and  the  result- 
ant signal  is  correlated  with  another  pulse  of  the  same  dur- 
ation. The  resultant  operation  is 

- J rft(i)  ■*  { orc*y  PunciienJ  m (. TT  + jt ) <±k  (93) 


which  in  the  frequency  domain  is 


c(f)  • W)£e*  M (P) 


A- < 


Jtr(f + C0»  9 


(94) 
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j is  simply  the  transform  of  the  autocorrelation 
function  of  the  pulse.  Thus, 


W)|J  - TW(f T)  5 T 


(95) 


When  inverse  transformed  this  function  is  a triangle  of 
height  T and  width  2T  (Fig.  A-l).  This  then  is  considered 
the  ideal  undistorted  signal  waveform  to  be  passed  through 
the  array  for  detection.  Eqs  ( 9^)  and  (95)  were  used  as  the 
basis  for  the  computer  program  that  was  written  to  investi- 
gate and  demonstrate  the  properties  of  the  array. 

Computer  Program  Analysis 

The  program  was  written  to  be  quite  versatile  in  its 
applications  for  this  study.  A series  of  specialized  sub- 
routines were  written  to  perfortmaost  of  the  work.  The  main 
program,  THESIS,  was  basically  a series  of  calls  to  these 
subroutines  to  obtain  a variety  of  results  for  visual  anal- 
ysis. The  major  subroutines  used  were  SIGNALS,  ARRAY,  and 
FFT.  The  supporting  subroutines  were  FREQ,  GRAPH,  ROTATE, 
DFT,  and  IDFT.  The  Fast  Fourier  Transform  (FFT)  algorithm 

used  relies  on  the  data  list  being  some  positive  integer 
M 

power  of  two,  2 . Thus,  all  the  subroutines  were  written  to 
conform  to  this  specification  for  vector  sizes. 

Subroutine  SIGNALS  fills  a vector,  Sigspoc,  with  equally 
spaced  frequency  samples  of  the  right  side  of  Eq  (95). 

Since  the  frequency  spectrum  of  Eq  (95)  is  infinitely  wide, 
the  frequency  spacing  must  be  chosen  so  that  an  adequate 
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Frequency  Domain 


amount  of  the  spectrum  content  is  sampled  to  yield  accept- 
able signal  waveshapes  in  the  time  domain  after  inverse 

p 

transforming.  This  particular  waveshape  (sine  x)  has  been 
studied  extensively  and  it  is  well  known  that  sampling  al- 
gorithms that  include  samples  from  one  or  two  of  the  side- 
lobes  yield  adequate  results  for  most  operations.  For  this 
program,  the  frequency  spectrum  is  sampled  so  that  samples 
from  at  least  the  seventh  sidelobe  on  either  side  of  the 
main  lobe  are  included  in  Sigspec. 

Subroutine  ARRAY  implements  the  complex  summation  in 
Eq  (9*0.  However,  the  form  of  the  summation  was  altered 
slightly  before  implementation.  An  f was  factored  from 
(f  + fQ)  leaving  fQ(l  + f/fQ)  . Next  r was  redefined 
as  r - dl  Aa  where  A,  is  the  wavelength  of  fQ.  Thus, 

d,  is  the  distance  r;  measured  in  units  of  A.  . Substi- 

•Hi  + iryr  Tnooo  ol.an^Juu  n +■  i ^r|  (CM  ol  rio 

£ geos*  _ 

*•1  Jxi 


■l  e^(' 


) d;  cos  O 


(96) 


because  c . Finally,  fQ  was  chosen  arbitrarily  to 

be  350  MHz,  a UHF  frequency.  Subroutine  ARRAY  then  fills  a 
vector,  ANTENNA,  with  equally  spaced  frequency  samples  of 
Eq  (96).  These  are  the  same  frequency  samples  used  in  sub- 
routine SIGNALS. 

The  third  major  subroutine  is  FFT.  It  is  a Discrete 
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Fourier  Transform  (DF?)  that  uses  the  decimation-in-time 
principle  (Ref  13 « 290-302)  for  computing  the  DPT  of  vectors 
that  are  a power  of  two  in  element  length.  The  routine  was 
written  so  that  either  a discrete  transform  or  a discrete 
inverse  transform  could  be  obtained  by  calling  FFT. 

These  three  subroutines  then  became  the  workhorses  for 
the  computer  analysis.  SIGNALS  was  called  to  create  the 
signal  spectrumt  ARRAY  was  called  to  create  the  antenna 
spectrum,  and  a simple  term  by  term  product  of  the  two  re- 
sulted in  the  array  output  vector  (Fig.  A-2).  This  output 
vector  was  then  passed  to  the  FFT  for  inverse  transforming. 
The  resultant  signal  was  then  graphed  so  that  it  could  be 
visually  compared  with  the  ideal  output  signal  of  Figure  A-l. 

The  five  supporting  subroutines  were  written  to  support 
the  details  of  graphing  the  output  of  the  three  main  subrou- 
tines. FREQ  has  the  simple  function  of  placing  the  appro- 
.prHa+.ft  values  in  the  horizontal  axis  vectors.  It  creates 
the  appropriate  time  and  frequency  vectors~l;hat  match  the 
sampling  values  used  by  SIGNALS  and  ARRAY.  GRAPH  has  the 
function  of  graphing  the  output  in  an  acceptable  format  by 
using  basic  CALCOMP  plotting  routims.  This  routine  is  only 
useful  with  the  facilities  available  at  FIT'S  computer 
center. 

ROTATE  is  a subroutine  that  was  used  in  the  analysis  of 
the  linear,  equally  spaced  array.  Section  III.  C.  proved 
that  the  signal  and  its  first  time  derivative  were  90#  out 
of  phase.  ROTATE  determines  the  phase  of  the  signal  and  then 
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shifts  the  phase  of  the  output  so  that  the  signal  can  be 
found  on  the  real  axis.  The  first  time  derivative  is  then 
located  on  the  imaginary  axis.  In  this  format  these  two 
quadrature  signals  can  easily  be  graphed  separately. 

DPT  and  IDFT  were  written  to  support  FFT.  DFT  and  IDFT 
set  the  parameters  of  FFT  to  cause  FFT  to  calculate  either 
the  discrete  transform  or  the  discrete  inverse  transform, 
respectively. 

All  the  above  subroutines  were  designed  so  that  they 
need  no  changes  for  any  run.  Obviously,  subroutine  SIGNALS 
would  have  to  be  rewritten  if  a different  input  signal  is 
desired. 

The  main  program,  THESIS,  is  written  v/ith  all  the  para- 
meters for  various  situations  included.  The  program  was 
written  for  batch  processing  and  little  consideration  was 
given  to  optimizing  this  routine.  The  various  parameter 
cards  were  simply  replaced  with  new  cards  that  contained  the 
new  parameters.  A listing  of  the  most  often  used  program  is 
included  as  Figure  A-3* 

Discussion  of  Data 

Any  time  a digital  computer  is  used  to  simulate  a con- 
tinuous system  errors  are  introduced.  These  errors  must  be 
identified  and  accounted  for.  This  computer  program  intro- 
duces errors  into  the  output  from  three  major  sources.  All 
three  are  discussed  here. 

The  first  effect  is  roundoff  error.  The  computer  per- 
forms calculations  to  a certain  degree  of  accuracy.  Rarely 
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Figure  A-3.  Fortran  Program. 
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Figure  A-3.  Fortran  Program. 
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Figure  A-3.  Fortran  Program. 
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i 8 this  precisely  the  right  number.  Under  moBt  circumstan- 
ces this  causes  no  real  problems.  However,  if  the  correct 
answer  is  in  itself  a very  small  number,  then  the  error  can 
be  large  compared  to  the  actual  answer.  This  effect  was  ob- 
served in  the  Tables  of  Chapter  III  where  A and  B were  tab- 
ulated. The  computer  never  calculated  their  values  to  be 
precisely  zero  at  locations  where  it  was  shown  that  they  are 
zero. 

Roundoff  errors  caused  other  problems)  problems  that 

affected  the  output  of  the  graphics  routine.  To  illustrate 

the  biggest  effects  of  roundoff  we  will  use  Figure  28  arid 

Figure  43.  Figure  28  shows  that  the  output  is  zero.  The 

program  did  not  calculate  all  the  output  samples  to  be  pre- 

cisely  zero.  They  had  values  that  ranged  from  10  to 
~28 

-10“  . When  the  graphics  routine  went  to  graph  these  points 

it  scaled  the  axes  so  that  the  details  of  the  "noise"  were 
visual.  We  had  three  possible  ways  to  deal  with  this  situa- 
tion.  (1)  Let  the  system  graph  the  "noise".  After  a few 
runs  we  decided  this  method  was  visually  unacceptable. 

(2)  Artificially  set  the  range  of  the  scaling  so  that  the 
"noise"  is  graphed  as  zero.  This  was  a possible  route. 

(3)  Pick  a value  that  is  considered  the  noise  threshold. 

Any  values  below  this  will  be  replaced  by  zero.  We  took 
this  option  because  it  was  so  simple  to  implement.  However, 
it  must  be  emphasized  that  we  only  performed  this  replace- 
ment at  the  time  of  graphing.  The  calculated  values  were 
used  as  is  for  all  operations.  Figure  43  shows  the  other 
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effect  that  the  roundoff  errors  caused.  The  roundoff  noise 
causes  the  phase  to  be  calculated  as  slightly  less  than  90° 
or  slightly  greater  than  90°  . This  combined  with  the  sec- 
ond source  of  computer  error  to  produce  the  results  of  Fig- 
ure 43.  The  arctangent  function  interprets  results  only  in 
the  range  from  90 * to  -90°  . Thus,  a value  slightly  greater 
than  90°  is  viewed  as  a value  slightly  less  than  -90*  . 

Thus,  the  noise  around  90*  or  -90*  of  phase  shift  causes 
wild  fluctuations  in  the  graph  of  Figure  43.  There  were 
several  methods  thought  up  to  artificially  eliminate  this 
graphics  problem,  but  they  were  not  implemented. 

The  third  error  deals  with  the  discrete  nature  of  the 
FFT.  The  FFT  is  a truncated  Fourier  Series  representation 
of  the  data.  In  order  to  include  enough  frequency  harmon- 
ics to  make  the  square  wave  look  exactly  square  would  re- 
quire a large  value  of  FMAX  and  the  number  of  frequency 
samples  2n.  This  is  due  to  Gibbs  phenomenon  which  states 
that  signals  with  discontinuities  cannot  be  represented  by 
finite  length  Fourier  Series.  The  result  of  a finite  Fourier 
Series  representation  is  the  oscillation  near  the  discon- 
tinuity that  appear  in  Figures  A-4,  A-5»  and  A-6  (Ref  2»105- 
107).  Figures  A-4,  A-5»  and  A-6  illustrate  how  the  oscil- 
lations on  the  square  wave  decrease  as  the  frequency  range 
and  the  number  of  samples  is  increased.  Figure  A-4  used  a 
256  point  FFT  and  only  sampled  frequency  components  up  to 
2 MHz.  Figure  A -5  used  a 512  point  FFT  and  sampled  to  4 MHz. 
Finally,  a 1024  point  FFT  and  a maximum  frequency  of  8 MHz 


130 


were  used  to  obtain  Figure  A-6.  These  figures  clearly  illus- 
trate that  the  higher  harmonics  aid  in  making  the  square  wave 
more  square.  The  plots  in  the  main  text  all  used  the  1024 
point  FFT. 

Summary 

This  appendix  has  presented  the  details  of  the  computer 
program  that  was  used  to  obtain  most  of  the  plots  presented 
in  the  text  of  this  thesis.  It  has  described  the  equations 
that  were  programmed  and  the  actual  Fortran  program  has  been 
included.  As  a final  step  a discussion  was  given  of  the 
effects  that  the  digital  program  had  on  the  form  of  the 
output.  This  included  the  observed  effects  of  roundoff 
errors,  the  limitations  of  the  arctangent  routine  used,  and 
the  observable  effects  of  using  a finite  length  Fourier 
Series  representation  of  waveforms. 
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